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Abstract. We rigorously prove existence and conformal invariance of scaling limits of magnetization and multi-point spin 
correlations in the critical Ising model on an arbitrary simply connected planar domain. This solves a number of conjectures 
coming from physical and mathematical literatures. The proof is based on convergence results for discrete holomorphic 
spinor observables. 



(N- 

^ ■ 1. Introduction 

. The Ising model plays a central role in equilibrium statistical mechanics, being a standard example of an order-disorder 
D ■ phase transition in dimensions two and above. Besides purely mathematical interest, its studies have found successful 

, applications to ferromagnetism, lattice gases, chemical adsorbates, ecology and image processing. 
fT^ ■ The phase transition in the Ising model in two dimensions has been a subject of extensive study, both in mathematics 
] and physics literature. Onsager [Ons44] exactly computed the free energy and the critical exponents of the model. Further, 

■ a gradual understanding of the Ising model at criticality led to the conjecture by Belavin, Polyakov and Zamolodchikov 
(~| ] that its scaling limit (as well as scaling limits of other critical models) should be conformally invariant and described by 
Q_( Conformal Field Theory [BPZ84a, BPZ84b]. 

In the mathematical language, this conjecture may be formulated as follows. Consider a continuous planar domain 
and let fls be discrete approximations to this domain. Then, as the mesh size 6 tends to zero, some quantities (expectations, 
probabilities etc.) in the model, under a proper normalization, have conformally invariant or covariant limits. Moreover, 
G ■ Conformal Field Theory predicts exact formulae for these limits. 

[ The mathematical proofs of such predictions were out of reach until recently. In 2006, Smirnov [Smi06] has rigorously 
^-H ■ established conformal invariance of the fermionic observables in the Ising model. Later, this result has been proven to be 
K*" I universal [ChSml2], and led to the proof of convergence of the interfaces in the Ising model to Schramm's SLE3 curves 
00 ■ [CDHKSll]. The full-plane energy correlation functions (that is, the correlations of n pairs of neighboring spins) were 
I computed on periodic isoradial graphs by Boutillier and de Tiliere [BoDT09, B0DTO8], and the scaling limit for bounded 

■ domains has been rigorously treated in [HoSmlO, HonlOa], confirming the CFT predictions for the energy field. 
The corresponding question about the spin correlations remained open, and the existing results mainly concern the 

case of the full plane. The diagonal and horizontal spin-spin correlations in the plane were computed by Wu [McWu73]. 
' Recently, Pinson [Pinll] has proven the rotational invariance of the full-plane two-point correlations. A number of 
remarkable results were obtained for the massive limits, see [WMTB76, SMJ80, PaTr83] and references therein. Palmer 
[Pal07] justified the CFT predictions at criticality by taking the zero-mass limit. In the presence of boundary, however, 
even the simplest problem of the asymptotics of magnetization of the Ising model in the half-plane remained open since 
it was posed by McCoy and Wu [McWu67]. 

In this paper, we rigorously prove existence and conformal covariance of scaling limits of all multi-point spin correlation 
functions in any simply-connected planar domain with + boundary conditions. For example, our result for magnetization, 
when combined with the asymptotics of the two-point correlations in the plane [McWu73], reads as follows. Let E^_^ [aa] 
denote the magnetization (expectation of spin) at the cite a with + boundary conditions in a discrete domains fls- Then 

6-^E+^[aa]'^Ci-Ta.d-i{a,n), 

as fig approximates f2 and the mesh size S tends to zero. Here Ci is an explicit lattice-dependent constant, and rad(a, fl) 
denotes the conformal radius of fl as seen from a, or, in other words, rad(a, 51) — \ip'(0)\, where 93 is a conformal map 
from the disc {2; G C : |z| < 1} to 51 mapping the origin to a. 

We establish similar results for the multi-point correlations K'^ [caoCai • ■ •"'ofc], and, in the case of two-point function, 
for free boundary conditions (see Theorems 1.1 and 1.3). Our previous results [Chlzll] immediately allow one to treat 
alternating +/— boundary conditions as well (see Remark 1.4). The method we use also applies to mixed correlations. 
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involving spin, energy, disorder and boundary change operators, and extends to multiply-connected domains, which will 
be worked out in a subsequent paper. 

The explicit formulae for the the scaling limits of E'^^[aao(Tai ...,cro^] were derived by the Conformal Field Theory 
methods in [BPZ84a, Car84, BuGu93]. In our approach the answers come in an explicit, but different form, being defined 
via solutions to a special interpolation problem, so one should additionally check that they coincide with the known CFT 
predictions. In this paper, we do this check for k = 1,2,3 and for the case when all the points are on the same vertical line 
in the half-plane (see further discussion in Sections 2.7 and the Appendix), which gives all the properties of the continuous 
correlations we need in the proof. 

Our main tool is discrete holomorphic fermionic observables. We use the spinor version of those, introduced in [Chlzll]. 
Fermionic observables per se were essentially invented back in 70's [KaCe71], and reminiscents of the spinor versions can 
be found in the works of Kyoto's school [SMJ77, SMJ79a, SMJ79b, SMJ80]. However, rigorous proofs of convergence 
require the delicate analysis of discrete boundary value problems for the observables, developed more recently [Smi06, 
SmilO, ChSml2]. 

Simultaneously and independently of our work, Dubcdat derived the 2n-point spin correlation functions in the full plane 
at criticality, combining exact bosonization techniques [Dubllb] and results on monomer correlations in the dimer model 
[Dublla]. Analogous results have been announced for polygonal domains. 



1.1. Main results. The Ising model on a graph ^ is a random assignment of ±1 spins to the vertices V of Q. The 
probability of a spin configuration a G {±1}^ is proportional to 

e-/3H(cT)^ 

where ^ > is the inverse temperature and H (a) := X^^^y cr^Ty is the energy of the configuration. We will consider the 
Ising model in discrete planar domains i^g which are subsets of the refining square grids rotated by 45° of diagonal mesh 
sizes 26 (thus, the distance between adjacent spins is V^S, see Figure 1). In this paper we are interested in the setup when 
ris approximate some simply connected domain 51 as 5 — )■ 0. 

If we consider the Ising model on the faces of fig with -|- boundary conditions (the spins on the boundary of il are set 
to -|-), there is a nattiral mapping (known as the low-temperature expansion, see [Pal07, Chapter 1]) between the spin 
configurations and the collections of closed contours (interfaces) on Jl^: trace an edge between any two faces with different 
spins. Under this mapping, the probability of a collection of interfaces lo C fls becomes proportional to 

Q,#cdgcs(w) 

where a = exp(— 2/3). The same procedure applied to free boundary conditions (no restrictions are set for boundary 
spins) produces collections of contours which also may start and end on the boundary of Cl. 

From now on we only consider the model at its critical point, which for the square grid corresponds to the parameter 
value /3c = ^ In (v^ + l), or equivalently to ac = — 1, when the phase transition occurs. As a normalizing factor, we 
will use the spin correlation in the unique infinite-volume limit of the critical planar Ising model: set 

q{5) :=EcJcro5(7iJ, 

where denotes the square grid S{l+i)Z'^, while 0^ and 1^ stand for proper approximations of the points 0, 1 e C. We 
first state the convergence theorem for the two-points functions, both with -|- and free boundary conditions: 

Theorem 1.1. Let V, be a bounded simply connected domain and Qg &e discretizations ofCl by the square grids of diagonal 

mesh size 2d. Then, for any e > 0, we have 



[^g^i.] ^ , , l^'^^b] ^ , .free 

— '-^^ > {(^aCTb)n and — '—^ > {(^a(Tb)^ ■. 



as 5 ^ Q, uniformly over all a,b G fl at distance at least e from dfl and from each other, where {cra(Jb)^ and {(TacrbyQ^ 
are explicit functions given by (1.1) and (1.4)- 

Remark 1.2. (i) Here and below we follow the traditional notation developed by physicists and denote the limits of spin 
correlations by the symbols {cFaa<^ai ■ ■ ■<^ak)Q, which should be understood as functions of the points ao,ai, . . . ,ak G 
which may also depend on O and boundary conditions b as parameters. 

(ii) The functions {aaa-b)^ and (cTaO'b)^**** are conformal covariants of degree | with respect to a and b: for any conformal 
map ip : ^ Q', one has 

(CTaab)f^ = {a^(a)<:r^{b))n' ■ \(p' {a)\^\(p' {b)\k (1.1) 

(iii) A celebrated result of T. T. Wu [McWu73] gives the rate of decay of g{d): 

where the lattice dependent constant is given by C2 = 2ie~^'''^~^\ and denotes the derivative of Riemann's zeta 
function. However, we do not need this result in our proofs. 
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We have a similar result for multi-point correlations with + boundary conditions: 

Theorem 1.3. Let SI be a bounded simply connected domain and Qs &e discretizations of SI by the square grids of diagonal 
mesh size 26. Then, for any e > and any k = 0,1, . . . , we have 

as 6 ^ 0, uniformly over all ao,ai, . . . ,ak € O at distance at least e from dSl and from each other, and the functions 
{(^ao^ai ■ ■ ■ <^ak)a havc the following covariance under conformal mappings (p : fi — > fi'; 

{cTao'^ai •••CTaJn = (o'y(ao)'^y(ai) " • " ' U.j=o\v' {dj) \^ ■ (1-2) 

The convergence of (fc + l)-point correlation with free boundary conditions (which by symmetry is only non-zero for 
odd k) can be obtained by our techniques as well, but for simplicity we do not include this in the present paper. 

Due to conformal covariance (1.1), (1-2), it is sufficient to give the explicit expressions for the correlation functions in 
the upper half-plane H. For magnetization and two-point functions, they read 

K)h = .J^ ,1 . (1-3) 



^"""^^^ = i2^a)H2^mb)r = i2^raa)H2^rab)i ' 

A different way to write these formulae involves conformal radius rad(a,f2) and hyperbolic distance d^: 

= 2ivadia,nyK 

(craCTfe)j, = ■ —Y, {(^a(^b)Q " 



(1.4) 



(l-cxp(-2d«(a,fe)))^' " exp(id«(a,fe))' 

For the discussion of the explicit form of ((Jq,, . . . (T,,^)^'^, fc > 2, we refer the reader to the Appendix. 

Remark 1.4 {alternating boundary conditions). Using results of [Chlzll], one can immediately generalize Theorem 1.3 to 
the alternating +/— boundary conditions. Namely, let b = . . . ,b2n} be a collection of points on dSl appearing in 
counterclockwise order and let bs C dSls denote a proper discretization of the family b. Let Eq^ denote an expectation 
for the Ising model (as above, defined on the faces of Slg) with the + boundary conditions on the counterclockwise arcs 
[b2j-i,b2j] and — boundary condition on the complementary arcs [&2j. ^'2j+i], j = l,...,n (we set 62n+i '■= bi). It 
was proven in [Chlzll] that all the ratios EQ^[o-ao. • ■^ak]/^ns{'^ao- ■ -^k] have explicit conformally invariant scaling limits 

Pj2(ao, . . ., Ofc). Thus, Theorem 1.3 implies that the expectations Q{S)~^EQ^[aaQ- ■ -CTak] themselves have conformally 
covariant scaling limits 



1.2. Key steps in the proof. In this section we list the key results that allow us to prove Theorems 1.1 and 1.3. We 

start our analysis with the discrete logarithmic derivatives of the spin-spin correlations with + boundary conditions. The 
following Theorem, which identifies the ratios of those correlation for neighboring faces is a cornerstone for the whole 
paper (keep in mind that our square lattice is rotated by 45°): 

Theorem 1.5. Let SI be a bounded simply connected domain and Sis be discretizations of SI by the square grids 5{l+i)l?. 
Then, for any e > and any k = 0,1, . . . , we have 

" 1 \ 

-11 '^tAn{a;ai,...,ak), (1.5) 

-1 -^mAn{a;ai,...,ak) (1.6) 




as (5 — > 0, uniformly over all faces a,ai, . . .Uk € Slg at distance at least e from dSl and from each other. The function 
An{a; ai, . . . , a^) is defined explicitly via solution to some special interpolation problem (see further details in Section 2.5) 
and has the following covariance property under conformal mappings ip -.Sl^ SI' : 

An{a;ai,...,ak) = An'{'p{a);p>{a{), . . . ,ip{ak)) ■ <p'{a) + ^^7^- (1-7) 

Proof. The proof is based on the convergence results for the discrete spinor observables. A rather delicate analysis is needed 
since we are interested in the values of observables near their singular points. See further details in Sections 2 and 3. □ 
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Now one can reconstruct the quantities log{aa<Tai ■ ■ ■ <^ak)ti integrating their derivatives ^n(a; ai, . . . , Ofe) inside f2. 
In particular, the constants in the integration can be chosen so that the exphcit formulae (1.3) and (1.4) given above hold 
true. In Section 2.8 we also discuss how to choose those multiplicative constants coherently for all and k in order to 
give a consistent definition of the continuous correlation functions. Note that the conformal covariance degree | in (1.2) 
is a direct consequence of the covariance rule (1.7) for vAn(a; 04, ... , aj.) (see Remark 2.21). 

Theorem 1.5 allows us to prove the following weaker form of the convergence result for the spin correlations: 

Corollary 1.6. Under conditions of Theorem 1.5, for any k > there exist some normalizing factors gk^i[S,Qs) that 
might depend on fl and the number of points a,ai, . . . ,ak but not on their positions such that 

as 6 ^ 0, uniformly over all faces a, ai, . . . afe e at distance at least e from dSl and from each other. 

Proof See Section 2.8. □ 

We now focus on the special case k = 1, i.e. on the two-point functions. The next theorem is a crucial tool which 
allows us to compare the normalizing factors g2{S,fls) with the full-plane case. We denote by Eq^^ the expectation for 
the critical Ising model defined on the vertices of (with free boundary conditions). 

Theorem 1.7. Let Q be a bounded simply connected domain and be discretizations ofCl by the square grids 6{l+i)Z'^. 
Then, for any e > 0, we have 



Bn{a;h), 



uniformly over all faces a,b at distance at least e from dil, and from each other, where ,8a (a; 6) is a conformal invariant 
of {^,a;b) which can be explicitly written via two-point functions (I.4) as Bn{a]b) = {aaO-b)^'^^ / {o-aO-b)^. 

Proof. The proof is based on the convergence results for the discrete spinor observables and the Kramers- Wannier duality, 
see further details in Sections 2.4, 2.6 and Section 3. □ 



The FKG inequality gives EqS'^ [cTo+aafc+a] < Ec^ [(^aC^b] < [ffaCfc] and it is easy to check that Ba{a;b) — )■ 1 as 6 
approaches to a. Since the normalizing factors g2{S,^s) do not depend on the positions of a, 6 e Q, we conclude that 
Eq^ [aaCTf,] ~ Ecj [cafTb] in the double limit when i5 — >■ and (later on) b a, thus relating g2{S,fls) with the full-plane 
normalization and establishing Theorem 1.1 (see further details in Section 2.9). 

Having the universal two-point function normalization g2{S,fl5) ~ g{5) as (5 — >■ 0, we derive asymptotics of all other 
gk+i{S, rii) using the following observation: as the point a approaches dfl, the continuous correlation functions are defined 
so that one has 

{(^a(^ai ■ ■ ■ CTaJn ~ (c^a)nKi • • ■ C^a Jo (1-8) 

and the same decorrelation result E^^ [o'aO'ai ■ ■ ■ o'ak] ^ E^^ [(Ta]E^^ [caj . . . aa^] holds true in the double limit 5 — )■ and 
a — >■ (see details in Section 2.10). This implies the recurrent formula gk+i{S, ^s) ~ ^?l(^, ^s)Qk{S, ^s) for A; = 1, 2, . . . 
and, further, gk+i{S,fls) ^ g{5)^''~^^^^'^ for all k, thus establishing Theorem 1.3. 



1.3. Organization of the paper. The rest of this paper is devoted to the proof of Theorems 1.1 and 1.3. Section 2 
outlines the proof and contains all the main ideas. Details, especially those involving hard discrete complex analysis 
techniques, are mostly postponed to Section 3. The reader not interested in these details may restrict himself to Section 2 
only. 

We fix the notation in Section 2.1. The main tool of this paper, the discrete holomorphic spinor observables, is 
introduced and discussed in Sections 2.2 and 2.3. In Section 2.4, we prove that the ratios of spin correlations that 
appear in Theorems 1.5 and 1.7 can be expressed in terms of these observables. In Sections 2.5 and 2.6, we discuss the 
continuous counterparts of the discrete observables and state the convergence Theorems 2.15, 2.17 and 2.19, that easily 
imply Theorems 1.5 and 1.7. We derive the formulae (1.3) and (1.4) in Section 2.7, and prove Corollary 1.6 in Section 2.8. 
We complete the proofs of Theorems 1.1 and 1.3 in Sections 2.9 and 2.10, respectively. 

Section 3 is devoted to the proof of Theorems 2.15, 2.17 and 2.19. We discuss the discrete properties of our observables 
and their full-plane analogue in Sections 3.1-3.3, and finish the proof of the main convergence theorems in Sections 3.4 
and 3.5. In the Appendix, we explain how to compute explicitly the continuous spinor observables in the half-plane. 
We also check that the correlation functions coincide with the CFT predictions at least when all points are on the same 
vertical line in the upper half-plane, and prove decorrelation identity (1.8). 
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Figure 1 . An example of a discrete domain ilg and notation for the sets of vertices {Vq^ ) , faces (Vq_^ ) , 
edge midpoints (Vq^) and (four types) of corners (V^_^ = V^^ U Vq^ U V^^ U Vq^). The mesh size i5 is a 
half-diagonal of a square face, thus the distance between adjacent spins is y/25. The inner vertices, faces, 
edges and corners are colored black, while the boundary ones are colored gray. 

Acknowledgements. All three authors wish to thank Stanislav Smirnov for introducing them to the Ising model and dis- 
crete complex analysis, for sharing his ideas and for constant support. We also would like to thank Stephane Benoist, David 
Cimasoni, Julien Dubedat, Hugo Duminil-Copin, Pierluigi Falco, Christophe Garban, Barry McCoy, Nikolai Makarov, 
Charles Newman, John Palmer, Yuval Peres, Mikhail Sodin, Tom Spencer and Yvan Velenik for many helpful discussions 
and comments. 

The work of Dmitry Chelkak and Konstantin Izyurov was supported by the Chebyshev Laboratory (Department of 
Mathematics and Mechanics, St. -Petersburg State University) under RF Government grant 11.G34.31.0026. Dmitry 
Chelkak was partly supported by P.Deligne's 2004 Balzan prize in Mathematics. Clement Hongler was supported by the 
National Science Foundation under grant DMS-1106588 and by the Minerva Foundation. Konstantin Izyurov was partially 
supported by ERG AG CONFRA and the Swiss National Science Foundation. 

2. HOLOMORPHIC SpINORS AND CORRELATION FUNCTIONS 

2.1. Notations. We start by fixing the notation which is used throughout the paper. 

2.1.1. Graph notation. Recall that we work on the square grid rotated by 45° 

Cs '■= {S (1 + i) (m + in) : m,n € Z} . 

The mesh size S is hence the size of a half-diagonal of a square face. We often will identify the vertices of with the 
corresponding complex numbers, the faces of Cs with their centers, the edges of <Cs with their midpoints, etc. 

We call discrete domain of mesh size S a union of grid faces (see also Figure 1 for notation given below) . We say that 
ris is simply connected if the corresponding polygonal domain is simply connected. 

• We denote by IntV^^^ the set of all faces belonging to ils (which are identified with their centers), by IntV^^^ the 
set of all vertices incident to these faces, and by IntV^^^^ the set of all edges incident to Vq^ , which are identified 
with their midpoints (or the vertices of a medial lattice). 

In order to simplify the presentation, we also assume that fig has no fiords of a single face width, i.e., all the edges 
joining vertices from V^^ belong to . This technical assumption can be easily relaxed, if necessary. 

• We denote by dV^^ , dV'^^ and dV^^ the sets of boundary faces, vertices and edges, i.e. those faces, vertices and 
edges which are incident but do not belong to IntVj°^ , IntVj' and IntV^ , respectively (see Figure 1). 

. We set V^^ := IntV^^ U SV^, etc. 
Below we also need to work with four corners of a given square face separately (see the crucial Definition 2.1 below). 
For a given vertex v £ Cg, we identify the nearby corners with the points v ± and w ± ijz on the complex plane. 

• We denote by the set of all corners incident to the vertices from IntV^^^. We also set Vfj™ := V^^. U V™^. 
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2.1.2. Double covers. In this paper we often deal with holomorphic functions (both discrete and continuous) which are 
defined on a double cover of a planar domain f2 and changes the sign between sheets (i.e., have the multiplicative 
monodromy —1 around branching points). We call such functions holomorphic spinors. The following notation will 
be used below: 

• For a planar domain J7 and a G we denote by [fi, a] the double cover of fi \ {a} branching around a. All such 
double covers are naturally viewed as subdomains of [C, a]. We often identify points on a double cover with their 
base points (so each base point is identified with the two points above it). 

• For several marked points ai, . . . ,ak G fl, we denote by [O, a; Oi, . . . , a/j] the double cover of f2 \ {a, ai, . . . , Uk} 
which branches around each of the points a, ai, . . . , aj,. 

• We will often compare spinors defined on [f2, a; ai, . . . , a^] with those defined on [C,a] (e.g., with 1/y/z — a or 
^/z — a) near the common branching point a. We write such equations meaning that they are valid in a small 
e-noigliborhood of a, with the natural correspondence of sheets. 

We will also consider double covers of discrete domains. In this case the following slightly modified notation will be 
convenient: 

• For a discrete domain fls and a face a G Int V^^ , we set [fls,a^] : = [Qs > \ { a + f } : excluding both points over the 
corner a+| from the natural double cover branching at a. Similarly, we set a~^] {oj}] := [^5, «; + |}) 
if several additional faces ai, . . . ,ak G IntV^^ are marked. 

2.1.3. Contours. In order to define the main tool of the paper - discrete holomorphic spinors, we need some additional 
notation related to the contour representation of the Ising model: 

• We denote by Cq^ the family of all collections of closed contours on fJ^, i.e. the family of subsets of edges ui C Vq^ 
such that every vertex v G Vq^ belongs to an even number of edges in w. 

Recall that Cn^ is in a natural 1-to-l correspondence with the spin configurations of the Ising model on the faces of fig 

with + boundary conditions (all the spins corresponding to boundary faces dV^j^ to +): trace an edge between 

any two faces with different spins. Below we also introduce the families of contour collections which, besides a number of 
closed loops, contain a single path running from one fixed corner x to another corner or an edge midpoint y: 

• For x,y € Vq™, let TTx,y = xvi . . . VnV be some walk with wi G Vq^ . We set Cn^ (x, y) := {w ® 7r,;,y, uj G Cn^ } , 
where © denotes the XOR, or symmetric difference. It is easy to sec that C^,, (x. y) does not depend on the 
particular choice of iix,y Note that, for any 7 G C^^ {x,y), there exists a (non-unique) decomposition of 7 into a 
collection of disjoint, simple loops and a path p (7) C 7 running from x to y. By a decomposition we mean that 
each edge in 7 G Cq^ {x, y) belongs to exactly one loop (or to p{'j)) and is visited only once, and that there are no 
transversal intersections or self-intersections (see Figure 2). 

2.2. Construction of discrete spinor observables. Now we are ready to introduce discrete spinor observables. The 
following definition generalizes the construction given in [Chlzll] to the case when a "source" point is inside Sis- 

Definition 2.1. Let Clg be a discrete domain and a,ai, . . . ,ak € IntVQ^ be inner faces. For a corner z G Vjq^ a^-ai ak] 
(below we also extend this definition to edge midpoints, see Remark 2.2(iii)), we define 



Fp„a;a,,...,a,] {z) ■= . ^ T E af'^^^^'^^ ■ <i>a;a.,...,a, (7, z) , (2.1) 



where 



#edges (7) is the number of full edges contained in 7 and ac = — 1 as above; 
the complex phase (l)a;ai,...,ak{l-:Z) is defined by (see also Figure 2) 

(7,^) := e-twind(p(7)) . (_i)#ioops„,a. '^^^^^^ • sheet„;„„...,„, (p(7),^), 

where, for a decomposition of 7 mentioned in Section 2.1.3, 

— wind (p (7)) is the total winding of the path p (7) when going from a + | to that is, the number of left 
tiurns minus the number of right turns, multiplied by |. 

— #loops„.„^_ (7 \ P (7)) is the number of loops in 7 \ p (7) that contain an odd number of marked points 
a,a\, . . . ,ak inside (equivalently, that change the sheet of the double-cover [fi^, a; ai, . . . , ttfe]), 

— the last factor sheeta;ai,...,afc (p (7) , z) is equal to +1 if z is on the same sheet of [fl^, a; ai, . . . , afc] as the lift 
of the end of p(7), and to —1 otherwise (more precisely, we fix one of the two points lying over the "source" 
a 4- 1 once forever and identify all other z G [Vis, a-., ai, • • • , cifc] with paths running from this a-|- | to 2; modulo 
homotopy and an appropriate index 2 subgroup of the fundamental group) ; 

the normalizing factor Z'^^ [ca^ai • • •O'afe] is defined by 

Z+ [a„<7„,...a„J := a^^^-^-) (-1)#'°°p^— -^('^^ 
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Figure 2. A contour collection 7 £ Cn^ (a + |, z) decomposed into non- intersecting loops (dashed) and 
a path p(7). Running from a + | to the projection of z, this path makes a Stt turn counterclockwise, thus 
g_±wind(p(7)) ^ rpj^g^g ^ gj^gjg j^^p ^ surrounding a, hence (-l)i°°P>'a(7\p(7)) = _i. Being lifted 
to the double cover [51^, a], this path ends on the other sheet, thus sheeta(p(7), z) = —1, and ^0(7, z) = —i. 



Remark 2.2. (i) For any collection of closed contours lu £ Cq^, the sign (-l)^'°°P'»^"i coincides with the product 

of spins aaCai ■ ■ ■<^ak ill the corresponding Ising model configuration with + boundary values. Since q,^°'^sos(w) ^^^^ ^j^^ 
Ising weight of u, one concludes that 

where Z^^ = X^weCn af'^'^^''^ is the partition function of the model. In particular, Zj|^ [(ya<7ai ■ ■ ■ ^'a^] > 0- 

(ii) It is easy to check that the complex phase (f>a;ai,...,ak ili z) is independent of the choice of the decomposition of 7 into 
a path ^(7) and a collection of loops, e.g., see discussion in [Chlzll]. Note that there are four types of corners: lying to 
the right of a nearby vertex v, below v, to the left of v, and upper v. For each of these groups, the total turning of the 
path p(7) is defined uniquely modulo 2tt. Therefore, the discrete spinors introduced above always have purely real values 
at the first group corners, are collinear to A :— e^* for the second group, etc. This motivates the following notation: 

• we partition the set of all corners into four subsets Vq^, Vq^, V^^ and Vq^ depending on the position of a 
nearby vertex v g V'^^ (to the left, upper, to the right, below) with respect to the corner. 

(iii) We extend Definition 2.1 to edge midpoints z G Vj^^ a^-ai a^] adding the factor (cos^)^"'^ to the formula (2.1), 
with #edges (7) being the number of full edges contained in 7 and the complex phase 4>a;ai....,akilT z) being defined as 
above. Note that in this case each edge midpoint z can be reached by a path p(7) from two opposite sides. Hence, the sum 
in (2.1) contains both real and imaginary terms, thus producing a complex number as the spinor value -F'[ni.a:ai,....afc] (z). 

2.3. S-hoIomorphicity and boundary conditions. A version of discrete holomorphicity, the notion of s-holomorphicity 
was introduced in [Smi06] together with the nonbranching version of discrete holomorphic observables, as a tool to study 
the critical Ising model on the square lattice. The properties of such functions were further investigated in [ChSml2] for a 
more general class of graphs. On the square grid, s-holomorphic functions may be thought of as (more classical) discrete 
holomorphic functions whose real part is defined on Vq^ and imaginary part on Vq^ , extended in a particular way to 

Vq^, Vq^, and further to V™^ (see more details in Section 3.1). Our definitions resemble those in [Smi06]. 

Definition 2.3. With each corner a; G (with r G {l,i,X, A}), we associate the line £ {x) := tM in the complex plane, 
and denote by P((^x) projection onto that line, defined by 

We say that a function F : V^™ — > C is s-holomorphic in ils if for every x G Vq^ and z G V[™ that are adjacent, one has 

i^(x) = P,(,) [F{z)]. 

For functions defined on double covers, we introduce the notion of s-holomorphicity exactly in the same manner. 
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The following proposition contains the crucial properties of the observables -F[na,a;ai,...,afc] that will allow us to analyze 
their scaling limits. For z G 9Vq^, let fout(-2) denote the "outer normal to the boundary at z": the edge whose midpoint 
is z, oriented towards the exterior of the domain and viewed as a complex number. 

Proposition 2.4. The function -F[n5,o;ai,...,aji] is s-holomorphic and has (multiplicative) monodromy —1 around each of 
the marked points a,ai, . . . ,ak, thus being a discrete s-holomorphic spinor on a^-ai au] " 

[F[n„a;a^^...,a,]{z) V^^out(^)] = for all z e dVf^^^,^.,^_,^y (2.2) 
Proof. We give a proof (which is mainly based on the standard XOR bijection, cf. [ChSml2]) in Section 3.1. □ 

Remark 2.5. The boundary conditions (2.2) are a priori not robust enough to pass to the scaling limit: even if the limiting 
domain 17 has a smooth boundary, the discrete normal I'out {z) can possibly admit only the values \/2e*T and v^e^"^, 
and so does not (pointwise) converge to its continuous counterpart. These conditions become much nicer, if one finds a 
way "to integrate" the square of f[nj,a;ai,....afc] : the real part of this integral satisfies Dirichlet boundary conditions on dO. 
due to (2.2). This approach is not as straightforward, since the square of a discrete holomorphic function is, in general, not 
discrete holomorphic, and so does not have a well-defined discrete primitive. However, it was noted in [Smi06] that one 
can naturally define the real part of the integral, using the s-holomorphicity of discrete observables, which is a stronger 
property than the usual discrete holomorphicity. Moreover, a technique developed in [ChSml2] allows one to treat this 
real part essentially as if it were a harmonic function, see further details in Section 3.3. 

2.4. From discrete spinors to ratios of correlations. The following lemma expresses the ratios of spin correlations 
(or their discrete logarithmic derivatives) from Theorem 1.5 in terms of the spinor observables introduced in Section 2.2. 
Moreover, it also provides the crucial ingredient for the proof of Theorem 1.7. 

Lemma 2.6. For any A; = 0, 1, ... , we have 

En, [<ya+2S(Ta,...(Ta^] („,3S\ 

Ta.+ r r- - ^[nii,a;ai,...,ak][0'+ -2) , (^•<3) 

where we take the comer a + ^ on the same sheet as the "source" point a + |. Moreover, for k = 1, we also have 

En, [f^ao-bj 

where E^S** denotes the expectation for the critical Ising model defined on the vertices of Q.s ( with free boundary conditions 

5 

outside the set Vq^) and the sign ± depends on the sheet where the corner b + ^ is taken. 
Proof. Recall that 

E#odges(7) , / I 3i5\ 

p I _ 7eC.,(a+f.a+f)«c (7, ^ + ^) 

^[Qs,a;ai...,ak] v'^ + 2 / ~ 'r+ r 1 ' 

while, taking into account Remark 2.2(i), 

E^^ [aa+2SO-ar ■ • ■ Q-g J ^ Z^u-£C», (_U 

There exists an obvious bijection between the sets Cn, (a + |, a + ^) and Cn,: removing the two corner-edges (a, a + |) 
and {a + 5,a + ^) from a given contour 7, we obtain a contour a;(7) € Cn, and vice versa. So, it suffices to show that 

0a;a....,a. (t, « + f) = (-l)#^°°P«<-+--i . 

Let us pick any loop in 7 and remove it. The left-hand side (respectively, the right-hand side) has changed the sign if and 
only if there was an odd number of points a,ai, . . . .a^ (respectively, a + 26, ai, . . . , a^) inside the loop. However, no loop 
in 7 separates a form a + 26 (such a loop would intersect ^(7)), so the two sides can only change sign simultaneously. Thus 
it is sufficient to consider the case when 7 is just a single non-self-intersecting path p(7) running from a -|- | to o -|- 
which is treated by the following observations: sheeta;ai....afc (p(7), a + ^) = —1 if and only if there is an odd number of 
points a, ai, . . . .a^ inside the loop ^(7), and wind(p(7)) = 2n mod 4n if an only if uj('y) separates a from a + 26. 
For (2.4), the Kramers- Wannier duality (e.g., see [Pal07, Chapter 1]) implies 



En, [o-ao-b] 2^, [(TaO-b] 

hence it is sufficient to prove that the (purely imaginary) number (j)a;b {l,b + ^) does not depend on 7. We have 
#loopSa.j(7 \ p(7)) = 0, since any loop in 7 either surrounds both a,b or none of them (otherwise it would intersect 
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the path p(7) joining a + | and b + |). Further, let tt^^ = bvi . . . VmO> be some simple walk with vi,. . . ,Vm G V^^ . Then, 
^(7) U (6 + |, 6) U 7r°j, U (a, a + |) is a loop, which we denote by /(7). Let n{-y) be the number of self-intersections of ^(7) 
(i.e., the number of intersections of p(7) with 7r^„). Since exp[— |wind(Z(7))] = (— 1)"'(t)+i (this is true for any closed 
loop), we see that 

exp[-iwind(p(7))] = • exp[iwind((6 + |, 6) U U (a, a + f))], 

where the second factor docs not depend on 7. But wc may also view tt^^ as a cut defining a sheet of the double cover 
[Q, a; b], meaning that sheeta;b(]3(7), + |) = (— 1)"^^-', which proves the desired result. □ 

Remark 2.7. Similarly to (2.3), one can check that 



= e^^ F^ns,a;au...,a,]{a+{iH)S)- (2-5) 



The proof boils down to the identity 

e±*<^a;a„...,a. (7, a + (l± f) j) = (-1)#'°°P^<'™— . 

for the natural bijection 7 i-)- ^(7) removing the two corner-edges (a, a + |) and {a + 5,a+ (1±^)(5) from a given contour 
7 € (a -|- |, a -|- (1 ± 5)^), which we leave to the reader. 

2.5. Continuous spinors. In this Section, we introduce the continuous counterparts of the discrete spinor observables 
defined in Section 2.2: the continuous holomorphic spinors f[Q,a;ai,...,ak]- We define them as solutions to the conformally 
covariant Riemann boundary value problem (2.6)-(2.8), which is a continuous analogue of the corresponding discrete 
boundary value problem (see Remark 2.10 below). 

Definition 2.8. Let be a bounded simply-connected domain with smooth bomidary, and a, ai,...,afe S fi. We 
define f[n,a;ai,...,ak] to be the (unique) holomorphic spinor on a; oi, . . . , Uk], branching around each of a, ai, . . . , and 
satisfying the following conditions: 



Sm [/[0,a;ai,...,a,](2)Vt'out(z)J = 0, Z G dfl; (2.6) 

Um sjz-a^ ■ fin 

,a;ai ,...,afcj 

(z) e iM.; (2.7) 



Im - a ■ fii2Miai,...Mk]iz) = 1- (2.8) 
where i'out(2) denotes the outer normal to the boundary of at z. 

Remark 2.9. (i) Note that a solution to the boundary value problem (2.6)-(2.8) is unique, if exists. Indeed, if /i, /2 denote 
two different solutions, then the spinor /i — /2 satisfies two first conditions, while lim^-j.^ \/z^^ ■ {fi{z) — f2{z)) = 0. 
Applying the Cauchy residue theorem to the single-valued function {fi{z) — f2{z))'^, one arrives at 

0<i-^ [ ihiz) - h{z)fdz = 2^ V hm {z - aj){h{z) - f2{z)f < 0, (2.9) 

where the first inequality easily follows from (2.6) and the second from (2.8). 

(ii) It is easy to see that if ly? : — ^ fi' is a conformal mapping, then we have 

1/2 

,a;ai . . . ,afc ](z) — I\ip{U),ip{a);ip(ai)...,ip{a,,)\ 

Mz))-{^'{z)r\ (2.10) 

We use this covariance property to define the continuous spinor in arbitrary simply connected domain 17, using a conformal 
map to some smooth bounded SI'. 

(iii) In Section 2.7 and the Appendix we give explicit solutions to (2.6)-(2.8) in the upper half-plane, thus proving its 

existence. 

Remark 2.10. The first condition (2.6) in Definition 2.8 is a natural counterpart of (2.2). The second condition (2.7) 
comes from the following observation: a discrete primitive -ff[o,o;ai,...,ofc] of the "differential form" 5Re [fj^^ ^.^^ ak]^^] 
(which may be defined due to the s-holomorphicity property of the discrete observable, see Remark 2.5) remains bounded 
from below near the branching points ai, . . . , Cfc as 5 — >■ 0. Thus, we impose the same condition for the scaling limits, 
which means that 3?e [/ f^^ ^.^^ ak]^^^ should behave like Cj log \z — aj\ for some negative cj as z — )■ aj, implying (2.7). 
The third condition (2.8) which fixes the behavior of the continuous spinor near the "source" point a is the most delicate 
one and will be clarified later on (see Lemma 2.13 and Section 3.2, particularly Lemma 3.5). Note that it is sufficient to 
assume that /[n,a;ai,...,ofc] does not blow up faster than l/y/z — a at a. Indeed, in this case the argument similar to (2.9) 
shows that lim y/z — a ■ fm a-ai ak]{^) > the rest is just the proper choice of normalization. 

Let us now introduce the quantities that play a central role in the computation of scaling limits of spin-spin correlations, 
appearing as the limits of their discrete logarithmic derivatives in Theorem 1.5. 
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We define the complex number An{(i', cti,. ■ ■ , ak) 

1 



Definition 2.11 

expansions 

f[n ,o;ai,...,aje 

near tlie point a. For the special case k 



A[n,a;ai,...,ak] ^s the Coefficient in the asymptotic 
+ 2^[a,„;„„...,„,]N/?^ + 0{{z - a)3/2) (2.11) 



\/z- a 

1, we also define the quantity BQ{a;b) 
near b: 



^^[Q,a;b] > as the coefficient in the 



vz — b 



(2.12) 



asymptotic expansion of the spinor f[n,a,ai,...,ak 

/[a,a;6] 

where the sign ± depends on the sheet of [O, a; 6]. 

Remark 2.12. Note that the covariance rule (1.7) for >l[n,o;ai,...,ofc] directly follows from the conformal covariance (2.10) 
of spinor observables: if <^ : f2 ^ f2' is a conformal mapping, then one has 



/[0,a;ai, 



,ip{a);ip{ai),...,ip(ak)] 



{z-a){l + §-^Az-a) + ...) 



1 



1 + 2 L4<^ • (a) + 



1 ^" (a) 



if{z) - <^(a) 
[l + 2A^-^p'{a){z 

{z-a) + ... 



[l + 2A^-{>p{z)-ip{a)) + ..:\ 
a) + ...] 



8 (a) 

where A^p = A[qi ^^(a);ip(a^),...,<p{ak)]- Similar arguments show that the coefficient B\Q^a;h\ is conformally invariant, i.e., 

Bn{a;b)=Bn'{'fi{a);ip{b)). (2.13) 

Thus, it is sufficient to find those quantities for some canonical domain, e.g., for the upper half-plane H. This is done in 

Section 2.7 for fc = 0, 1 and in the Appendix for fc > 2. 

2.6. Convergence of spinors. The main purpose of this section is to give a rough sketch of the proofs of Theorem 1.5 
and Theorem 1.7. The complete proofs are given in Section 3 below, which is the most technical part of our paper. We 

use the following conventions concerning convergence of discrete s-holomorphic functions; 

• we say that a family of discrete domains (fia)^ approximates a continuous domain C C as 5 — >■ 0, if converges 
to dO, in the Hausdorff sense (note that our proofs can be easily generalized for the Caratheodory convergence of 
planar domains which is weaker than the Hausdorff convergence used in this paper for simplicity); 

• we say that an s-holomorphic functions (or a spinors) defined in fi^™ (or its double cover) tend to a holomorphic 

function (or a spinor) / as 5 — )• 0, if the "mid-edge values" Fs I approximate the values of /, while the "corner 

' "j 

values" F\ , t G {1, A, i, A}, tend to the projections of / onto the corresponding lines tM (see Definition 2.3); 

• we say that a convergence of discrete functions Fs{z) = Fs{z; a, 6, ... ) to f{z) — /(z; a, fe, . . . ) is uniform on some 
compact set, iff the differences \Fs{z\ a,b, . . .) — f{z; a,b,...)\ are uniformly small as 5 — > 0, when we interpret 
lattice vertices (or mid-edges, corners, etc) z,a,b,... as the corresponding complex points when we plug them 
into /. 

The crucial ingredient of our proofs is the interplay between (a) the values of discrete spinor observables near their 
branching points, which are related to the ratios of spin correlations by Lemma 2.6, and (b) the mid-range behavior of 
these observables, which can be further related to the asymptotics expansions of their scaling limits (2.11), (2.12). 

As a main tool to relate (a) and (b), we use a full-plane version F^c^ a] of the spinor observable (since we are interested 
in local considerations, it is sufficient to stick to the case of one marked point). In principle, it can be constructed as an 
infinite-volume limit of the finite-domain observables, but we prefer a more descriptive strategy, which is outlined after 
the following lemma claiming the existence of -F[C5,o]- 

Lemma 2.13. For a G Vq, there exists a (unique) s-holomorphic spinor F^c^^a] '■ ^[c"a->] ^ ^ ^^^'^ that F[c^^a]{0'+^) = 1 
and F{z) = o(l) as 2 —> oo. Moreover, 



1 



= : /[C,a](^)) 



(2.14) 
(2.15) 

discrete harmonic 

x,x < 0}. This definition is motivated by the following 
observation: in the continuous setup, the function Jfte Wj^/z — a\ can be viewed as the properly normalized harmonic 
measure of the (small neighborhood of) tip a in the slit plane C\{a-|-a;, x < 0}. Second, we extend -F[C5,a] to Vq by harmonic 



s->-o ^z — a 

uniformly on compact subsets ofC \ {a}, where the normalizing factor ■&{6) is defined as 

m :=F[c„„](a+f +2^L^J). 

Proof. The detailed proof is given in Section 3.2. First, we define the (real) values of F^c^.a] '^Cs 
measure of the tip point a -|- ^ in the slit discrete plane \{c 
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conjugation, then by symmetry to another sheet of [C5,a], and eventually as an s-holomorphic function to V^^,, and 
Vq. The convergence (2.14) follows from known results on convergence of the harmonic measures [ChSmll]. □ 

Remark 2.14. The normalizing factor 'd{S) is essentially the value of F^Cs,a] at a + 1. In Section 3.2 we show that 

C-\/S <^{5) <C+\/S, (2.16) 

where C± > are some absolute constants. Note that one can compute the limit lim5_>.o 'd{S)/\/5 using the recent work 
of Dubedat [Dub 11a], but we do not need this sharp result. 

We further use the normalizing factors 'd{6) introduced above in order to formulate the following convergence theorem 
for discrete spinor observables away from dfl and a, ai, . . . , a^: 

Theorem 2.15. Let discrete simply connected domains Clg approximate a bounded simply connected domain il as 6 ^ 0. 
Then, for any e > and any k = 0,1, . . . , we have 

';^J^F[Qs,a;ai,...,ak]{^) /[n,a;ai,...,afe] (•2), 

uniformly over all a,ai, . . . ,ak € Vq^ and z G Vq^ which are at the distance at least e from dCl and from each other. 

Proof See Section 3.4. □ 

Since we are interested in the coefficient A[Q^aiai,...,ak] in front of the term in (2.11), along with a discrete 

analogue -F'fCs.o] of the function I/1/2 — a given by Lemma 2.13, we also need a discrete counterpart Gfc^.a] of the function 
\J z — a. We construct Gjc^, by the "discrete integration" of -Fp^.a]! just mimicking the continuous setup. It is sufficient 
to define Gy^^^a\ on V^^ only, as Lemma 2.6 deals with the (real) values of discrete observables at the point a + ^ G V^.^. 

Lemma 2.16. For a G V^^, there exists a (unique) discrete harmonic spinor G[c^,a] • "'^[Q-a^] ~^ ^ vanishing on the 
half-line {a + x,x<0} such that G[c^,a]{z) = 0{\z — a\^/'^) as z^ 00, and G[Cs.a]{(i + ^) = Moreover, 

:^G[c,,a]{z) ^ Ucy/z-a =: g[c,a]{z) , (2.17) 

uniformly on compact subsets 0/ C \ {a}. 

Proof. For z G V^^ ^j, we set G[c^_a]{z) '■= ^ ' Sjlo F[Cs,a]{z — 2j(5) (certainly, one should check the convergence of this 
series and the harmonicity on {a + x, x > 0}). See further details in Section 3.2. □ 



In the continuum limit, the leading term in the expansion of f[Q.,a;ai,...,ak\'~ f[C,a\ near a is given by '2.A[Q„a;ai,...,akW z — a. 
It is hence plausible to believe that the same holds true for the discrete spinors, and one has 

{F[ns,a;ai,...,ak] " -^[Cj.a]) («■ + ^) ~ 25Re^[n5,a;ai,...,afe] " G[Cs,a\ (a + ^) 

up to higher-order terms (the real part appears due to discrete complex analysis subtleties, as real and imaginary parts 
of s-holomorphic functions are defined on different lattices, and a -|- ^ G ^n^)- justify this heuristics in 

Theorem 2.17. Under conditions of Theorem 2.15, we have 

F[n,,a;ar,...,a,\{a+^)-l-2'SlzAyn,a-a„...,a,]-5 = O {6) (2.18) 

as (5 — > 0, uniformly over all a,ai, . . . ,ak G Vq^ which are at the distance at least e from dfl and from each other. 

Proof. See Section 3.5. □ 

Remark 2.18. In the proof of Theorem 2.17 we show that F[na,a;ai,...,afe] and F'pa.a] are 5-close to each other at all points 
around a, in particular at a + (1 ± |)5. Together with (2.5), this implies 

^Oa [o'a-Kl±i)^0'ai • • • O'ofe] ^ ^ 

E^, [o-aO-oi • • • 0-aJ *^-0 

uniformly over all a,ai, . . . ,ak & Vq^ which are at the distance at least e from dCl and from each other. 

The similar analysis for the quantity B^d a-.h] which is defined, for fc = 1, by the asymptotic expansion (2.12), is even 
simpler since we need to match the first-order coefficients instead of the second-order ones. The result is given by 

Theorem 2.19. Under conditions of Theorem 2.15 for k ^ 1, we have 

F[n„a-M{b+l)±iB[a,a-M = o{l) (2.19) 
as S ^ 0, uniformly over all a,b G Vq^ which are at the distance at least e from dCl and from each other, where the sign 
± depends on the sheet of [ft, a; b] . 

Proof. See Section 3.5. □ 
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Proof of Theorems 1.5 and 1.7. Due to Lemma 2.6, asymptotics (1.5) is a reformulation of (2.18), while Theorem 1.7 
is equivalent to Theorem 2.19 (the sign ± is fixed due to the positivity of spin-spin correlations). To check (1.6), we rotate 
our domain O around a by 90° clockwise. According to conformal covariance rule (1.7), the coefiicient A multiplies by i, 
so the desired result follows from SRe [iA] = —QmA. □ 

2.7. Computations for one and two points. Now we illustrate the results obtained in the previous section constructing 
the continuous spinors solving the boundary value problem (2.6)-(2.8) for two simplest cases k = and k = 1. We use the 
upper half-plane H as a canonical domain. Since it is unbounded, one should impose the additional regularity condition 
at 00 ensuring that f[M,a;ai,...,ak] = 0{\z\~^) as ^ ^ oo, so that after a conformal mapping to a bounded domain with 
smooth boundary the spinor remains bounded. 

2.7.1. The case k = 0. For a single marked point a e H we have 

/[H.al(.)= ,^^'^""f_, , ze[C,a], (2.20) 
- a){z - a) 

since this spinor clearly satisfies both conditions (2.6) and (2.8). Writing down the asymptotic expansion (2.11) at a, one 
sees that 

Am{a) = -—z . 

Further, we define the quantity log(c7a)H as a primitive of the differential form ^e[AM{a)da] = -~dy/8y, if a = x + iy: 

(c7a)]j5j = const • (9m a)~4, a e H. 
The proper choice of the multiplicative constant which is discussed below leads us to the definition (aa)^ := 2 3 (23m a)~i . 

2.7.2. The case k = 1. Similarly, for two marked points a, 6 G H, one can check that the spinor 

(2z9ma)^ [{b -a){b- a)] Hz - b) + [{b - a){b - a)] i{z-b) ^ ^ 

/[H,o;6] (Z) = 77 — -I , |, — -T 7} =77 TT =rrrj^ , z e L, o; 6 , (2.21) 

\b-a\ + \b-a\ [{z - a){z - a){z - b){z - b)Y'^ 

satisfies (2.6)-(2.8). Looking at the asymptotic expansions (2.12) at b, one obtains 

(4Sma9m6)^ {agabyr 

BM[a; b) = jr — _, , r = -p. 

\b-a\ + \b-a\ {(Tacri,)^ 

Further, asymptotic expansions near a give 

A r u\- 1 \b-a\-\b-a\ 
- + 4{\b-a\ + \b-a\) 

Straightforward computations shows that one can define a primitive /" Jfe [>Ih(2^; b)dz] for a e H (here we treat b as 

a parameter) that coincides with the function log((TaCrb)g given by (1.4). Note that ((TaCrt,)^ is uniquely defined (up to 
a multiplicative constant) by its logarithmic derivative with respect to a and the symmetry {aaf^b)^ ~ (^b'^o)iH- T'^® 
multiplicative constants in (1.3), (1.4) are chosen so that {aacrh}^ ^ (cra)e (o"!))!! if) say, a — )• dM, resembling decorrelation 
properties of the critical Ising model which arc given in Section 2.10. Below wc show that, for all k, the similar integration 
procedure allows one to define real symmetric functions {o'aO'ai ■ ■ ■ o"afc)e such that 

(o-aCTai ...craj+ - (cra) + (cr„, ...craj + , a ^ dn. (2.22) 

2.8. Proof of Corollary 1.6 and decorrelation identities. Denote ao := a. The main purpose of this section is to 
define the functions {(JaoC^ai ■ ■ ■<^ak)n basing on their logarithmic derivatives with respect to each of aj and (eventual) 
decorrelation properties (2.22). Let us denote aj = Xj + iyj, and consider a differential form Cqj.+i on the manifold 
j^fe+i {(flQ, . . . ak) S 1^*^+^ : aj am,j < fn} of all {k + l)-tuples of distinct points in Q, defined as follows: 

k k 

-Ca.fe+i :=^SRe [AQ{aj;ai, ..,aj, ..,ak)daj] = '^[^eAQ{aj;ai,..,aj,..,ak)dxj-^mAniaj;ai,..,aj,..,ak)dyj] (2.23) 

J=0 ,;=0 

(we use the standard notation aj for the omitted argument). The following proposition is a straightforward corollary of 
Theorem 1.5 obtained via integration with respect to positions of points: 

Proposition 2.20. The form £o,fe+i is exact. If {aao ■ ■ ■ o'afc)n •= 6^P[/ '^n.fe+i] denotes the exponential of its primitive, 
then, as Sis approximates ft, one has 

En, Who 0-6i . . . CT6 J (o-fco (Tbi . . . 0-6 J ^ 



[<^ao • • • O-o J {aa^ O^^... ) ^ 

for all [k + lytuples {a^lg, e 0*^+^ 



(2.24) 
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Remark. Note that (2.24) is essentially equivalent to the claim of Corollary 1.6, and the proof also ensures that this 
convergence is uniform, if all aj are at distance at least e from dfl and each other, and the same holds true for bj. 

Proof. Color the faces of fls black and white, in a chessboard fashion. By Remark 2.18, the ratio of spin-spin correla- 
tions at two adjacent spins tends to 1, uniformly away from the boundary, so we can assume that all aj,bj are colored 
white. Let G Cl he such that [ao,aQ] is a horizontal segment contained in O and disjoint with ai, . . . ,afc. Denote by 
do = vo,vi, . . . , = a'o a straight horizontal lattice path approximating [ao, a'^]. Then, by Theorem 1.5, one has 

log— ^— r = — ^— r - 1 [l + o{l)) = 26- [mAa[vj;ai,...,ak) + o{l)\ 

as 5 ^ 0, where the o(l) terms are unform in j. Consequently, 

ka' CToi . . . (To J ^ r 

log— — -=2S-y JttAn{vj;ai,...,ak) + o{l) — > / ^eAa{x + iy;ai,...,ak)dx. 

A similar formula with — 3m^j2(a; + iy; ai, . . . , a^) in the right-hand side applies to the case when [ao, Qol ^ vertical 
segment. Moreover, one can move other points ai, . . . , at along horizontal and vertical segments as well. Hence, one has 

where 7 is any path in 17*^+^ that connects (A: -|- l)-tuples and and consists of segments (in 17*^+^) with all but 

one coordinates fixed. Since the left-hand side does not depend on the choice of the path, the form Cn^k+i is exact. □ 

Remark 2.21. Note that {cTaof^ai ■ ■ ■ (^ak)n is a symmetric function of ao, ai, . . . , afe since it is constructed as a primitive of 
the symmetric differential form jCa^k+i- Integrating the conformal covariance rule (1.7) for An, one immediately obtains 



logp2^^1111^!^ = f\z[An{z;a,,...,au)dz 

Waa^a-L ■ ■ ■ Cat/o Jan 



I 

J a, 



1 f""'" 

5ie [A^{Q){^p{z);^p{a{),...,ip{ak))^p'{z)dz\+- I 3fie 

^ Jan 



ao ^ J ao 

{(^v(bo)(^v{a^)■■■(^V(a^))^(n) , 1, \(f'{bo)\ 

log -I- - log ■ 



viz) . 



Iterating, one arrives at the following covariance rule for the ratios of correlations which is a weaker form of (1.2): 

{aboO-bi ■■■o-fcJn ^ (g'y(bo)^y(fci) ■■■g'y(bfc))y(n) ^ yr \v'{bj)\i ^5) 

{(TaoCTa, ■ ■ ■ C^ak)n ■ ■ ■ ^ v(o.k))'^{Q.) j=0 I («J ) I ^ 

Proposition 2.20 defines the continuous correlation functions {(Jao^ai ■■■^ak)n multiplicative constants, which 

may depend on Q. and fc, since the primitive of £n,fc-i-i is defined up to an additive constant. A natural way to choose 
these constants coherently for all domains and any number of points is suggested by the following lemma. Denote by 

the quantity that measures how deeply in the bulk of Q the points a,b are. Denote also {(Ja(Jb)c \o, — b\~^. 

Lemma 2.22. There exists a unique way to choose the multiplicative normalization of {cFaCFai ■■■<^ak)n so that for all 

domains ^l and for all k = 0,1, .. . the following holds: 

{cTaaa,...aak)^ ^ ^ (^a^fc)^ ^ ^ Da(a; 6) ^ 0. (2.27) 

(0"a)n(0-ai ...O-ajQ {(^a(^b)c 

Being normalized in this way, the correlation functions {aa(7ai ■ ■ ■<^ak)n satisfy the conformal covariance rule (1.2). 

Proof See Section A.3. □ 

Recall that the properties (2.27) are directly motivated by the corresponding properties of the discrete correlations, 
see Lemmas 2.24 and 2.26. Summarizing, we define the continuous correlation functions ((TaUai ...aak)n ^'^ 
exponentials of the primitives of £a,fe+i, normalized as in Lemma 2.22. As it was discussed in Section 2.7, in the 
particular cases = and = 1, this definition reproduces formulae (1.3) and (1.4), respectively. The case fc > 2 is 
discussed in the Appendix. 
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2.9. From ratios of correlations to Theorem 1.1. This section is devoted to tlie proof of Tlieorem 1.1. Our goal is 
to relate the normalizing factors P2{5, ^s) from Corollary 1.6 (which, in principle, might depend on f2), with the full-plane 
normalization p{5). The proof is based on Theorem 1.7, which claims the convergence of the ratios of free and + spin-spin 
correlations to an explicit limit Bn{a;b) = {cra<Jb)n^ / {o'a(Jb)^i- ^^^o classical FKG (e.g., see [Gri06, Chapter 2]) 
and GHS (see [GHS70]) inequalities for the Ising model. The small additional ingredient is given by 

Remark 2.23. The following is fulfilled: 

Bnia;b)^l as Dn(a;6)-)-0, and Ba(o;6)-)-0 as a ^ dQ, (2.28) 

where the quantity 'DQ{a-,b) is given by (2.26). This follows readily from the conformal invariance of Bn(a;6) and the 
explicit formulae (1.4) in the half-plane. Also, one has 

^^^^ ^ 1 as Vn{a; b) ^ 0. (2.29) 

Indeed, let <^ be a conformal map from f2 to HI such that (p{a) = i. By standard estimates for conformal maps, when 
X>n(a; 6) — )• 0, one has (p{b) — )• i and |<p'(6)| — |<^'(a)|. Hence, the explicit formula (1.4) imply 

^ {(^v,{a)'^vib))t¥{a)\~'\^'{b)\^ ^ \v{a)~^{b)\-i\ip'ia)\i ^ ^^^^^ ^ ^ 

(0-a0-&)c (o'aCTh)c \a-b\~i 

Lemma 2.24. For any t] > there exists an e > such that the following holds: ifDa^a; b) < e and Vts approximates ^, 
then 

l-r?<%l^<l 

provided that 6 is small enough. 

Proof. By FKG inequality, lE^j'^ [^aCfe] < lEc^ [caCfo] < [caCr;,] for any domain containing a, b, hence the right-hand 
side readily follows. For the left-hand side, choose A^ = $1* — 5. We have, by Theorem 1.7, 

K'rK'Tb] Ef5SqT„+,,T,+a] 

^ ^ ' ' Bn{a;b). 



E+ [aaab] E+ [a„(76] s^o 

Due to Remark 2.28, one can choose e so that Bn(a; 6) > 1 — ^, which gives the result. □ 

Proof of Theorem 1.1. Fix 77 > 0. For shortness, below we will write a + e for its lattice approximations a + 2^[^J, 
and a -f 1 for its lattice approximation a + 2S[jg\ . Note that 

Eg, [ggO-b] ^ E+^ [o-ggb] ^ [a gab] ^ [ggQ-g+e] ^ Eq KcTg+e] 

q{5) [o-gCTa+l] E^^ [cTaCTa+e] Ec^ [cTaCTa+c] Ecj [cTgCra+l] 

By Lemma 2.24, we can find a small e > and a large domain As containing a, b, a + e and a+1, such that 
1 . IE+ [aaaa+e] ^ ^Cs K^g+e] E+ [aaaa+i] 

— V ^ T — 1 — ?y < T ■ , < 1 -h ??, 

Eq [cTaCTo+eJ Eq [(Ta(Ta+iJ E^^ [aaO-a+e] 

provided that 6 is small enough. Consequently, 

_ .2 ^(is ^A, [o"g'^g+^] J^n, [g^g^b] ^ , n '^n, K^fc] [ffaCTa+e] 



^-Us [^a(Ta+r] El^ [cTgCTg+l] g{S) E^^ [craCTa+e] Ej^ [aaaa+l] ' 

and, by convergence of the ratios proven in Proposition 2.20, 

{aa(Ja+e)n (crgCrg+l)J ~ ffiS) ~ {aa(Ta+e)a (CTgCTg+l)A 

for 6 small enough. Since r] can be chosen arbitrary small, and the bounds do not depend on 6, it only remains to show 
that we can make the factor 

as close to one as wc wish by choosing e small enough and A large enough. However, this follows readily from (2.29) if we 
multiply this factor by {aa(Ta+i)c = 1- Thus, E'^^ [cTgCifc] ~ Q{5){aaCrb)Q as 5 -> 0. To derive the asymptotics of two-point 
correlations for free boundary conditions as 5 — >■ 0, note that Theorem 1.7 implies 

[<^a<Tb] ~ Bn{a; b) ■ E+ [cjaOb] ~ q{5) ■ Bn{a; b){c7aC7b)tt = Q{5){aa<JbYr ■ 

The fact that we have — 5 instead of 0,$ plays no role, since they both approximate the same continuous domain f2 and 
the convergence of {q{5))~^E'^ [(Ja(^b] is independent of the particular choice of lattice approximations. □ 
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Remark 2.25. As a simple byproduct of our analysis, we obtain the rotational invariance of the full-plane correlations 

recently proven by Pinson [Piull]: by FKG inequality, for any (large) domain fi^, one has E^^"" [o-aCb] < Ec^ [caCt] < 
^fis Wa'^b] and, due to Theorem 1.1 and (2.28), both sides have the same asymptotics when Qs exhausts Cs- Then, (2.29) 
gives the desired result: (p((5))^"'" Ec^ [cTaCfc] ~> {<^a<^b)c — W~ b\^i as 5 — > 0. 

2.10. Decorrelation near the boundary and the proof of Theorem 1.3. This section is devoted to the proof of 
Theorem 1.3. Note that it was already proven above in the special case fc = 1, as a part of Theorem 1.1. Our goal is to 
relate the normalizing factors in the Corollary 1.6 with g{5). Below we rely upon decorrelation identities (2.27). 

Lemma 2.26. Given a domain D, with marked points ai, . . . , Ofe and a number r] > 0, there is an e > such that the 
following holds: if a € fl is e-close to the boundary and Q,s approximates il, then 

^ E+ K]E+ K...a.J 



Proof. The upper bound follows readily from FKG inequality. For the lower one, consider first the case = 1. A celebrated 
application (see [DeMolO]) of the GHS inequality reads EQ^[crocr(jJ — E^^ [tTaJEQ^ [ctoJ < EQ^''[cracraJ, or equivalently, 

^ Efr-Ka.J ^ E+ HE+ K] 



As 5 — ^ 0, the left-hand side converges to 1 — BQ{a; ai), thus the claim follows from (2.28). 

To prove the result for A; > 2, assume that we have already proven Theorem 1.3 for all k' < k. Let 7 be a crosscut 
(simple path) in ^l separating a from ai,...,ak, and let r2(a) and f2(ai) be the corresponding connected components. 
Then E'^^^^s^ [a-a]E^^^i^^^^ [aai ■ ■ ■ (Tq^] is equal to the correlation of cto, CToi , • . • , ct^j. in fig conditioned on the event that all 
spins neighboring 7 are pluses. Hence, FKG inequality implies 

Eo,(a)K]En,(aoKi ■■■^g.] E+Jaa]E+ K ...g,,] _ E+^^^^^E+^^^^^K . . ^ ^ 

i'^aj^Q, [cTai • ■ ■ E+^ [aafJa, • • • J E+^ [aafJa^ • • • J 

By our induction assumption, the first term on the left-hand side converges to 

(^")n(a) (^ai ■ ■ ■ "'afc)n(ai) 
(O'ai)n {(^0.1 ■ ■ ■ (^ak)n 

hence it is sufficient to show that we can make this quantity arbitrary close to 1 by choosing appropriately a and 7. We 

first fix a crosscut 7 in such a way that 0(ai) would be Caratheodory close to ft as seen from ai and then put a deeply 
inside fla, so that fla would be Caratheodory close to f2 as seen from a. If two domains are Caratheodory close, then the 
conformal maps from these domains to H mapping the pole to i (and say, with positive derivative there) are uniformly 
close on compacts together with their derivatives. Thus, the lemma follows from continuity of the half-plane functions 
((Tai • • ■'^ak)u with respect to positions of ai, . . . , a^;. □ 

Proof of Theorem 1.3. Recall that the special case A; = 1 is already done as a part of Theorem 1.1. We proceed by 
induction which (together with the proof of Lemma 2.26 given above) starts as follows: Ti&Li Tq 1^2 ^2 is 
T^=> . . . (where Tk and Lk mean the particular cases of Theorem 1.3 and Lemma 2.26, respectively). 
Fix ?7 > 0. For any c G we can write 

/e+ _E+ [aJE+ [a,] E+ E+ [a,] 



Vv^/ E+K^,] q{6) E^ 

By Lemma 2.26, if we choose c close enough to the boundary, then 

l-r?< ^+ ^ <1, 

provided that 6 is small enough. Due to Theorem 1.1 and Proposition 2.20, one li; 



as 



Consequently, 
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provided that 5 is small enough. Since 77 is arbitrary, and we can make the term {cracrc)^/ {a)Q{c)Q arbitrary close to one 
(recall that this particular case of (2.27) follows from explicit computations in the half-plane given in Section 2.7), we 
finish the proof by remark that positivity of magnetization fixes the sign. 

To get the convergence of g{S)~^~E^[aaa-ai ■ c^at]) proceed by induction and write it as 

QiS)'^ KKcTa^ • • .^aj ' E+[(7e]E+[ag, . . . CTg J ' ^/^Q^iS) 

The proof is finished as before: take c close to the boundary, estimate the second term in the right-hand side by Lemma 2.26, 
then use convergence of the first and the third terms (which follows by Proposition 2.20 and induction hypothesis, 
respectively) and asymptotics (2.27) for the continuous correlations: {(7c<Jai ■ ■ ■ '^ak)n ^ {(^cj^i^ai ■ ■ ■ ^ak)tt as c — )• dO.. □ 

3. Proofs of the main convergence Theorems 2.15, 2.17 and 2.19 for spinors 

3.1. S-holomorphicity of discrete observables. The notion of s-holomorphicity was essentially introduced by Smirnov 
in [Smi06] and used for the study of the planar critical Ising model in [Smi06, SmilO, ChSml2, HonlOa, HoSmlO, Chlzll]. 
Our definitions resemble [HoSmlO] and are equivalent to those of [ChSml2, Chlzll] after the multiplication by see 
also Section 3.3.1 below. 

Recall that we associate the line i{x) = tM. in the complex plane with each corner x G (for r e {l,i,A, A}), and 
denote by Pi(^x) the projection onto that line: 

Pi{x) M = i (w + T^w) , yiv e c. 

Definition 2.3 says that a function F : Vq™ — ^ C is s-holomorphic in (we also use the same definitions for double 
covers) if for every x € Vq^ and z £ Vq^ that are adjacent, one has 

F(x) = P,(,) [F{z)]. 

Remark 3.1. The set V^;* = V;^^ U V^^ may be viewed as a square lattice, divided into V^. , V^. in a chessboard fashion. 
By definition, the restriction of an s-holomorphic function F to Vq* is real on and purely imaginary on . It is not 
difficult to check ([Smi06, ChSml2]) that this restriction is in fact discrete holomorphic, that is, for any a € Vq' one has 

i ■ [F{a +{1 + i)5) - F{a)] = F{a + id) - F{a + 6). (3.1) 
Moreover, the converse is also true: given a discrete holomorphic function F : V^'^ — >■ RU iR, we can first extend it to Vq 

by the formula F{x) := F{x + if) + F{x - j|) and then to V^f by F (x) := P^(^) [F [x ± |)] , since (3.1) ensures that 

the result is the same for both choices. 

We now check the s-holomorphicity of discrete spinor observables, essentially mimicking [ChSml2, HoSmlO, Chlzll]. 

Proof of Proposition 2.4- Let z G ^.^^ be a medial vertex and x e Vj^^ ^.^^ be one of four nearby 

corners (so that |a; — 2;| = |). We should check that 

Pi{x)[F[n^,a;ai,...,a^:]iz)] = F^Qs ,a;ai,...,ak] (x) , (3.2) 

where both sides are defined as sums over the contour sets Cq^ [a+ ^,z) and Cq^ (a -|- | , a;) , respectively. There is an 

obvious bijection 7 h- > lu{j) between these two sets provided by taking XOR (symmetric difference) of a configuration with 
two half edges (zv) and (vx), where v denotes the vertex which is adjacent to both x and z. Hence, it is sufficient to check 
that for any configuration 7 e Cq^ [a+ ^,z) one has 

(cos|)-i . P,(,) [af dges(^) . </.g,g„...,g,(7,z)] = af^"^'^)) • <f>a.a,....M7),x) (3.3) 

(the additional factor (cos f )~^ comes from our defintion of the discrete observable on medial vertices, see Remark 2.2). 
There are two cases: either (zv) is contained in 7, which leads us to 

#edges(a;(7)) = #edges(7) and exp[-|wind(p(a;(7)))] = 6=*=^ exp[-5wind(p(7))], 

or not, which gives 

#edges(a;(7)) = #edges(7) -|- 1 and exp[-|wind(p(a;(7)))] = e^^ exp[-|wind(p(7))]. 
Let us also note that in both cases 

(_1)#100P...„„...,„,(7) Sheetg,g„...,g,(7,z) = (-l)#'°°P^-l ».(-(^»ShCCtg^g,,...,g,(.;(7),z) 

since if oj destroys a loop that changed the sheet of lu (leading to the change of the first factor) , then this loop becomes a 
part of p(7), so the second factor changes simultaneously. Thus, one can factor out Q,#'''^ses(w(7)) _ (p^^^^^^^^^^^(^ijj(^'-f)^x) from 
both sides of (3.3). In the first case (3.3) readily follows, while in the second one it becomes equivalent to 

(cos f )"^ cos ^ = \/2 - 1 = ac 
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Thus, i^[o,a:ai ,...,afc] IS s-holoniorphic. The factor sheeta;ai,...,afc (7) z) in the definition guarantees that it has multipHcative 
monodromy ^1 around each of the marked points o, ai, . . . , Uk, and in order to prove that /[n,a,ai,...,afc] obeys the boundary 
conditions (2.2), it is sufficient to note that, if z is on the boundary, then wind(p(7)) = Voutiz) mod 27r. □ 

Note that the spinor F^^ a-.ai a^] is not s-holomorphic at a + |, as the identities for windings in the above proof fail. 
The next lemma shows, in particular, that its values at the nearby medial vertices a+ have different projections onto 
the imaginary line iM. = £[a + |), so one cannot extend i^[n,a;oi,...,ofe] to a + | in an s-holomorphic way. 

Lemma 3.2. For medial vertices a + taken on the same sheet of [fig, a;ai,..., ak] as the "source" a + |, one has 

PiR [F[n,,a;a„-,at] + ^^)] = =F«, 

Proof. We consider the point a+^-^, the other case is similar. Given a configuration 7 £ Cs{a + ^,a + 5^^) and applying, 
as above, the XOR bijection with two half-edges {a + ^,a + S) and {a + S,a+ ^^S), we obtain a configuration 0^(7) € €5. 
Since the normalizing factor Zq^ [a, ai, . . . , ak] is a sum over Cs, it is sufficient to show that, for any 7, 

(C0S|)-1 • P« [ap''^'''^^Ua;u^,...,udl,a+^S)] = -i(-l)#^°°P^-»i .,(-(7))„#edges(a,(7)) (3.4) 

Consider two cases, as in the proof above. If {a + 5,a+ ^^6) G 7 (respectively, {a + 6,a+ ^^5) ^ 7), then 

#edges(u;(7)) = #edges(7) (respectively, #edges(w(7)) = #edges(7) + 1) . 

Also, we may disregard all loops in 0^(7) that do not contain the edge {a + 6,a + iS), as they contribute the same sign 
to both sides. If {a + 6,a + ^-^6) € 7, then wind(p(7)) = ^ mod 47r, otherwise wind(p(7)) = ^ mod 47r. Hence the 
lemma boils down to the following elementary identities: 

(cos |)"i • PjR[e"^] = -i and (cos ■ P«[e"^] = -iuc- □ 

3.2. The full-plane discrete spinor F^Cs,a^] and its discrete primitive G[Cs,a'*]- This subsection is mainly devoted 
to the construction of the full-plane analogue of discrete spinor observables, as announced in Lemma 2.13. After this, we 
also prove Lemma 2.16 and the double-sided estimate (2.16) of the normalizing factor i?((5). 

3.2.1. Discrete harmonic measure in the slitted plane. We start with an important technical ingredient - the discrete 
Beurling estimate with optimal exponent ^. On the square lattice, it was obtained by Kesten [Kes87j, and then generalized 
by Lawler and Limic [LaLi04]. Let denote the slitted discrete plane VQ\ia5 where La := {a + x,x < |}. 

Lemma 3.3. For all 2; € and A c and some absolute constant C > 0, the following estimates are fulfilled: 

W<C<5^k-«r^ (3.5) 

hm|''(a+f) < C(55(dist(a;A))-^ (3.6) 

where hm^^ (z) denotes the discrete harmonic measure of the set A in viewed from z, i.e. the probability for the simple 
random walk started at z to reach A before it hits the boundary of Eg . 

Proof. This easily follows from [LaLi04] and simple reversibility arguments for random walks. □ 
Below we also need some additional estimates for the discrete harmonic measure hm^"* 3^ . and its "discrete derivatives". 

Lemma 3.4. (i) If z e Eg is such that^t{z — a) < and \'^m{z — a)\ < ||3?e(z — a)| (i.e., z is close to the cut La), then 

hmf ^,^{z)<CSi\'^m{z-a)\-\z-a\-^ . (3.7) 

(ii) For all edges {zz') G S5, one has 

<CS^z-a\-i. (3.8) 



(Hi) After the normalization by the value ■(?((5) at the proper lattice approximation of the point a + 1, the function 
■!?(5)~^hm|^^3iT. converges to 5ie [i-fs/z — a\ uniformly on compact subsets ofC\La. Moreover, this convergence holds 

true in C^-sense meaning that the "discrete derivatives" (3.8) converge to the corresponding partial derivatives as well. 

Proof, (i) To prove (3.7), note that the probability for the random walk started at z to leave the ball of radius — a| 

around z before hitting d'E^s is 0(|9m (z — a)| • |z — a|~^), and if that happened the probability to hit a+^■ is still uniformly 
bounded byO((5^|z — a|~^) due to the Beurling estimate (3.5). 

(ii) The estimate (3.8) for the discrete derivatives follows from (3.5), (3.7) and the (discrete) Harnack estimate (e.g., see 
Proposition 2.7 in [ChSmll]), applied to the ball of radius ^\z — a\ (or W^m {z — a)\, if z is close to La) around z. 
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(iii) This is essentially a special case of Theorem 3.13 in [ChSmll] which claims the convergence of discrete Poisson kernels 
normalized at some inner point to their continuous counterparts. The fact that our domain is unbounded plays no role here 
since the positive discrete harmonic function t?(5)~"'^hm|^^3a | is uniformly bounded in any annulus z : r < \z — a\ < 

< r < R. Indeed, if it is big at some point v in the annulus, then, by the maximum principle, it is also big along some 
nearest- neighbor path running from v to a + ^ . The discrete harmonic measure of such a path as seen from a + 1 is 
bounded from below by a constant, leading to a contradiction, since our function equals 1 at this point by definition of 
In particular, i?(^)^^hm|^_|_3j| is uniformly bounded near any circle around a, and so everywhere outside this circle 



too, thus preventing appearance of the term \/z — a in the limit (which corresponds to the mass at oo and plays a 
crucial role in the proof of Lemma 2.16 given below). □ 

3.2.2. Construction of the full-plane spinor fp^.a]- We are now ready to construct -F[C5,a]i as announced in Lemma 2.13. 

Proof of Lemma 2.13. Let S± denote two copies of Cs\La- We first define the real component F^^^ on Vl_^ as 

i^[ic„a]W=±hinp^^^(z), zeVh^. (3.9) 

Since F^^^ vanishes identically on La, by identifying the upper side of this cut in S+ with the lower side in S_ and vice 
versa, we obtain a function F^^^ which is defined and harmonic everywhere on except at the (two) points over 

a + ^. Recall that Lemma 3.4(iii) ensures the convergence 

which holds true uniformly on compact subsets of C \ Lai as well as the convergence of "discrete derivatives" of F^^ ^^{z) 
to the corresponding partial derivatives of 3?e [l/V^ — o]. 

We then define the imaginary component F^-^^ : V^^^ = V^^ a]\i'^ + |} — > iM as discrete harmonic conjugate of 
F^^^ ^j, that is, by integrating the identity (3.1) along paths on V^^ starting from, say, one of the two fibers of the 
point a + ^. Since is harmonic, its discrete harmonic conjugate is well-defined at least on the universal cover of 

[Ca, Further, let zu denote some simple loop in Vq, starting and ending at a + symmetric with respect to the 
horizontal line {a + x,x G R}, and surrounding the singularity (so, it lifts to [C5,a~*'] as a path connecting two different 
fibers of a + ^). It follows from the (anti)symmetry of F^^^ with respect to the cut La (note that zu changes the sheet 
once passing across La), that the total increment of along w is zero. Thus, vanishes at both fibers of the 

point a + and hence inherits the spinor property of its discrete harmonic conjugate F^^__ 

Note that F^^^^ vanishes not only at a + ^ which was chosen above as the starting point of the discrete integration, 
but also everywhere on the half-line Ra ■= {a + x,x > |} since, for 5 := a -|- | -|- 2j5, j > 1, the discrete holomorphicity 
equation (3.3) and the symmetry of with respect to the horizontal line {a + x,x gM.} imply 

F(cs,a]ib + 26) - F[^^_„](6) = z ■ [Fic„a](6 - iS) - 2F^Cs,a]ib) + Flcs,a]ib + (2 - t)6)] = AFic,,a](6 + ,5) = 0. 

Further, the uniform estimate (3.8) guarantees that F^^^ is uniformly bounded (just take a path of discrete integration 
in the definition of F^^^ running from a+ ^■ to \z\ and then along the circular arc to z). 

It is worthwhile to note that F^^^ also admits a discrete harmonic measure representation similar to (3.9). Namely, 
let Sj_ denote the two sheets of [C5,a~*^] \ Ra, where signs in the notation are chosen so that Fj^^ > in the upper 
half-plane on (correspondingly, in the lower half-plane on S'_). Then, 

Indeed, on each of the sheets Sj_ one can further extend F^^ (as a harmonic conjugate of F^^^ ) to the point a -|- | : 
the only obstruction to do so on [Ca, a] was that the increment of F^^^ along the smallest loop surrounding a+^- would 
be non-zero since FjJ,^ is not harmonic at a+^■, but on Sj_ this loop intersects the cut Ra- Therefore, the (bounded) 

function F^^^ is harmonic on Sj_ and vanishes on Ra\{a+^}, thus being proportional to \mi^^^^^[z). Finally, it follows 
from the symmetry arguments that 

hmj^^, J («+(§+ i) 5) - lim=^*+|j (a + |) = hm^^^^^ (« + + i) 5) - hm^^+^j (« + f ) 

= ^^[C„a](«+(5+^)'^)-^[C„a](«+¥) 
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which fixes the multiplicative constant. Again, Theorem 3.13 of [ChSmll] implies 

^/ica^) ^„ ^-7^^ (3-12) 

uniformly on compact subsets of C \ Ra, together with the discrete derivatives. Now. following Remark 3.1, we extend 
Pc5,a] to Vp^j as an s-holomorphic function. The convergence (2.14) readily follows from (3.10) and (3.12) (for points 
lying near La, one can approximate by summing discrete derivatives of and vice versa near Ra). □ 

The following lemma explains the future role of F^c^^a]'- it has the same "discrete singularity" (whatever it means) as 
F[ns,a,ai,...,ak] at a + |. This allows one to handle this singularity by taking the difference Fp^^a;ai,...,ak] - PlCs,a] 

(a) A 

Lemma 3.5. The spinor Fjj^^ ^ := F^Q^^a,ai,...,ak] ~ ^[Cj.a]) extended to be zero at a+ ^, is s-holomorphic in any 

simply- connected neighborhood of a which do not contain other branching points ai,...,ak- 

Proof. By Lemma 3.2, it suffices to show that P^r [Fp^^a] (a + "'r^'^)] = on the sheet S+. The considerations given in 
the previous proof show that -Fp^^o] can be extended in an s-holomorphic way at a + | such that 

P« [F[c,M [a + ^5)] = -i ■ hm j(a + f ) = -i, for a + ifi^ e 5'+. 

However, by definition, coincides with in the upper half-plane, but with S_ in the lower half-plane. □ 

3.2.3. Construction of the harmonic spinor G[Cs,a] ■ Now we "integrate" the real component of F[c^,a] in order to construct 
a discrete counterpart of the harmonic spinor 5f e \/z — a, as announced in Lemma 2.16. Along the way, we also prove 
the double-sided estimate (2.16) of the normalizing factor 'd{6). Note that the upper bound is given by discrete Beurling 
estimate (3.5), so we need to prove the lower bound only. 

Proof of Lemma 2.16 and the estimate (2.16). We use notation introduced above. Define, for z e V^, 

oo 

G[C„a](^) := 5 ■ E^[C.,a](^ - m- (3.13) 

(due to the estimate (3.7), this sum always converges). Clearly, G[Cf,.a\ vanishes on the left cut La- We are going to prove 
that it is harmonic everywhere inside S+ (including the point a-\-^■). For z outside Ra, the harmonicity of G[C5,o] readily 
follows from the harmonicity of F^^^ ^j, so let z e .Ro- Denote by Tliqiz) := V^^ H {z + w : —2N5 < 3?ew, Qmw < 2N5} 
the (sufficiently large) square centered at z and write the discrete Green formula: 

N 

E^^[Q,a](^-2j5)= E ^4.,alH= E (4.,al(^')-4.,alH)- (^.14) 

i=o weTiNiz) (ww')ednN{z) 

Let N be large enough so that |w — a| > N5 for all boundary edges {ww') of TZn{z). Then, it immediately follows from 
the estimate (3.8) that the last sum is 0{N ■ 62{NS)-i) = 0{N-2). Passing to a limit as TV ^ OO, we conclude that 

oo 

AG[c„a](^) = <5E AFic,,„,(^ - 2jS) = 0, 

so G[Cs,a] is indeed discrete harmonic in S+. Moreover, since it vanishes on the cut La one can harmonically extend G[Cs,a] 
to the whole double cover V^^^ by symmetry. 

Let ^{S) denote the value of G[C5, a] in a lattice approximation of the point a-|-l. Arguingasintheproof of Lemma 3.4(iii), 
we see that ^ 

—F¥:^lCs.a]{z) ^ 3fJe - a (3.15) 

(since G^c^^a] vanishes everywhere on La and remains bounded near by the maximum principle, in this case the limiting 
positive harmonic function should be proportional to 5Je ^/z — a, and the multiplicative normalization is fixed at a -\- 1). 
Moreover, this convergence is uniform on compact subsets of C \ io and the same convergence holds true for discrete 
derivatives. In particular, this implies 



thus yielding v{5) ^ as 5 — >■ 0. Note that this allows us to give a simple proof of the lower bound in (2.16): since 
discrete harmonic functions {v {S))^'^G[Cs,a\ are uniformly bounded near the unit circle around a and vanish identically 
on La, discrete Beurling estimate(3.6) implies 

W)- W) -""^^ 
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(and v{5) ~ '&{5) as i5 — )• 0). Finally, the convergence (3.15) near La follows from the convergence of a](-^)' since the 
tails in definition (3.13) are uniformly small due to (3.7). □ 

3.3. The boundary value problem for spinors. In this section we reformulate the Riemann-type boundary value 

problem for holomorphic spinors (both discrete and continuous) using primitives of their squares. Note that this approach 
is not completely straightforward, since the square of a holomorphic function, in general, does not have discrete primitive. 
However, it was noted in [Smi06] that one can naturally define the real part of the primitive, using the s-holomorphicity 
of observables. Moreover, a technique developed in [CliSml2] (see, in particular, sections 3.4 and 3.5 therein) allows one 
to treat this real part essentially if it was a harmonic function. Below, we summarize the tools we will use. We warn 
the reader that all our definitions are equivalent to those of [ChSml2] after the multiplication of the spinor by Vi, which 
means that imaginary part, sub-/super-harmonicity and positivity of functions and their (inner) normal derivatives used 
in [ChSml2] should be replaced by real part, super-/sub-harmomcity and negativity, respectively. 

3.3.1. Discrete integration of the squared spinor observables. Let A| be the standard (unnormalized) discrete Laplacian 
acting on functions Hg : Vq^ — > M (which are defined on faces of Sis)' 

A°gHl{z) := Y,{HI{w)-HI{z)), z G IntV^,, 

where the sum is over the four neighbors w E V^-^^ of z. Similarly, for functions H' : Vq^ — )• M defined on vertices of Q,s, 
let A* denote the slightly modified discrete Laplacian; 

A}H!{z) := ^c{w,z)iH'siw)-H!{z)), ^ € IntV^,, 

where the conductance c {w, z) is equal to 1 for inner edges (i.e., for w e IntV^^) and is equal to 2(\/2 — 1) for the boundary 
edges (see Section 3.6 in [ChSml2] or [DHNll] for the reason of this "boundary modification" of A*). 

Proposition 3.6. For an s-holomorph/ic spinor observable Fs = F^Q^^a;ai,...,ak] '■ [^<5: 'i^! ^^i, . . . , afc] — >■ C satisfying 
boundary conditions (2.2), one can define a real-valued function Hs = H[Qg a;ai,...,ak] '■ ~^ ^ which is a discrete 
analogue of the primitive 5fte / {Fs{z))^dz, so that the following properties are fulfilled: 

• for any adjacent w G Vq^ and v G Vq^ , one has 

H°s H - HI (v) = 25 \Fs {\{w + v))\^ , (3.16) 

where ^^{w + v) € Vq^^ is the corner between v and w (according to Lemma 3.2, we set \Fs (a + |) | := 1 ); 

• Hs satisfies Dirichlet boundary conditions: Hg (v) = for any w G QVq^ , and H* {v) = for any v G QVq^ ; 

• H* has a "negative inner normal derivative", i.e. Hs (u) < for any vertex v G V^^ adjacent to a boundary vertex; 

• Hg is Ag-subharmonic on V^j^ \ {'^^'^^ij • • • I'^k}, while H' is A'-superharmonic on Vq^ \ {a + 5}. 

Proof. All the claims follow directly from the results of Section 3.3 in [ChSml2j. Since all listed properties are local, the 
spinor nature of Fs plays no role here (note that the right-hand side of (3.16) does not depend on the sheet). □ 

Remark 3.7. By construction, Hg{w) > Hg{v) for adjacent w and v. Combined with super- and sub-harmonicity, this 
implies a maximum principle for Hs: if fl's C Cls does not contain a (respectively, any of a, ai, . . . , afe), then 

mm Hs = mm Hg (respectively, max Hs = max Hg). 

dQ', ^ Q', an',' 

Moreover, if h.mA{z) denotes the discrete harmonic measure of a set A in Q.'g viewed from z, then 

Hs{z)>{l- hm^(z)) minif* + \mY\{z) mini?*; 

cTi^ A 

Hs{z) < (1 — hm^(z)) ma.yiHg + hm^(z) mayiHg. 

Remark 3.8. The subharmonicity of Hg fails at a.a,i. . . . .Ok because i^[j2j,a;oi...afe] branches at those points, while the 
superharmonicity of H^ fails at a + 5 because of the discrete singularity of -F[n5,a;oi...afe] which is not defined at a + |. In 
particular, Lemma 3.5 implies that the function 

= ^[S,a;a,...,a.] ■= J {Fiiz))'dz, = ^^[n„a;a,...,a.] " %„a], 

(defined in the same way as Hs accordingly to Proposition 3.6) is superharmonic on Vq^ and subharmonic on V^^ \ {a} 
everywhere in a neighborhood of o. Moreover, since Fg^\a + |) = on both sheets, the values of Fg'^^ at the nearby 
corners a± |,a± ^ and midedges a± satisfy the s-holomorphicity conditions as if it were nonbranching at a. Thus, 

(a) 

Hg is subharmonic at the branching point a G Vq^ too. 
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3.3.2. Integration of squared spinors in the continuous setup. We now give a characterization of the continuous spinors 
solving boundary value problem (2.6)-(2.8) in terms of the primitives of their squares. In the next section, this charac- 
terization will be used in the proofs of main convergence results. 

Proposition 3.9. Let Q be a simply- connected domain, and suppose a holomorphic spinor f = .f[n,a:ai,...,ak] solves the 
boundary value problem (2.6)-(2.8) (or is defined according to Remark 2.9(ii), if Vl is not smooth). Define two harmonic 
functions 

/i := 3fJe/ {f{z)fdz and h^"^ := iRe / {f{z) - fcA^)fdz, 
where fc,a{z) ■= ^/y/z — a. Then, the following holds true: 

(1) h is a single-valued function in f2\{a, ai, . . . , Ok}, continuous up to dSl, which satisfies Dirichlet boundary condi- 
tions h = const on dO, (since h is defined up to an additive constant, we assume that h = on dCl); 

(2) there is no point zq on dfl such that h{z) >0 in a neighborhood of zq; 

(3) h is bounded from below in a neighborhood of each a\,. . . ,ak; 

(4) h'^^^ is single-valued and bounded in a neighborhood of a. 

Moreover, ifh and h^'^^ satisfy (l)-(4), then f solves the boundary value problem (2.6)-(2.8). 

Proof. Note that, being integrated, the covariance property (2.10) claims the conformal invariance of both h and h^^\ 
As the properties (l)-(4) are preserved under conformal mappings too, we will further assume that is smooth. Note 
that the property (2.6) is equivalent to (1) and (2): it states that f"^ {z)i'ont{z) is positive on the boundary, which yields 
^Taut{z)^ = (where t{z) denotes a tangent vector) and 9^„„t(z)ft- > everywhere on dfl. A straightforward integration of 
the asymptotics of / near aj and a given by the conditions (2.7) and (2.8), respectively, gives 

Kz) = -Cj\og\z-aj\ + 0{l), z^Gj, 
h^^\z) = log|0-a| +0(1), z^a, ^ ' 

for some Cj > 0, giving (3) and (4). Vice versa, (3), (4) imply (2. 7), (2. 8) by differentiating and taking the square root. □ 

3.4. Convergence of discrete observables away from singularities. In this section we prove the convergence of 
(properly normalized) discrete spinor observables to their continuous counterparts on compact subsets of fl\{a,ai, . . . , ak}- 

Proof of Theorem 2.15. Let the function 

2 



Hs := 5Re J (^.^Fs{z)^ dz, Fs = F[n„„;ai,., 



be defined on Vq° accordingly to Proposition 3.6. Given e > 0, define 

flsie) := flsn{z : dist(z; {a, ai, . . . , Ok}) > e}. 

Assume for a moment that, for any e > 0, the functions Hg remain uniformly bounded on ^^^(e) by some constant 
C(e) as 5 — >■ 0. Then by Theorem 3.12 in [ChSml2], the functions {'&{5))~^Fs are equicontinuous on ^^(e). Thus, by 
passing to a subsequence and applying the diagonal process, we can assume that (i!) {S))~^ Fs tends to a limit / and 
Hs h := J f^ uniformly on compact subsets of Sl\{a, ai, . . . , Uk}. Our goal is to check that / satisfies the properties 
(1) - (4) given in Proposition 3.9, then the uniqueness proven in Remark 2.9(i) imphes / = f[n,a;ai,...,ak]- 

Clearly, / is a holomorphic spinor on [O, a, ai, . . . , ak]- By superhamonicity of H', we have 

min Hs{z) = min H'{z), 

\z—aj\<2e e<\z — aj\<2e 

thus h is bounded from below in the neighborhoods of Oi, . . . , a^, so the property (3) holds true. By maximum principle 
for Hs (see Remark 3.7), taking into account that Hs = on dflg, we have that 

HUz) > -C(e)(l-hm^^(;)(^)) 
H°,{z) < C{e){l-hr4ll^\z)) 

Since hm^^^^^^(2;) — > 1 uniformly in ^ as 2; approaches the boundary of il^, this implies h = on dfl, giving (1). Moreover, 
thanks to Remark 6.3 in [ChSml2]. we also have (2): there is no point on dfl such that h > in its neighborhood. 
Consider now the following function, discussed (up to multiplicative normalization) in Remark 3.8: 

2 



ffi^) := 5Re I dz, = Fs - Fp,,] 



which is well-defined in the disc {z : \z — a\ < r} provided that r is small enough. Since both {{){6))^^Fs and {'d{6))^^F[Cij,a] 
uniformly converge to / and /[c,o] (see Lemma 2.13) on compact subsets of Aj. := {z : < \z — a\ < r}, Hg^^ converges 
to h^^^ := 5ie f^{f{w) — /[c,o](w^))^<^w; everywhere in Ar. By Remark 3.8 (sub- and super-harmonicity of Hg^^ on V° and 
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V* near the point a), the functions Hg^^ are uniformly bounded in Ar, so h^^^ is bounded in Ar too, which concludes the 
proof of the property (4). Therefore, / = f[n,a;ai,...Mk]- 

It remains to justify the uniform boundedness of Hg away from a,ai, . . . ,ak- On the contrary, suppose that 

Ms{e) := max lHsl — > oo 

for some e > and along some subsequence of ^'s. Then the (re-) normalized functions {Ms{e))~^Hs are uniformly bounded 
in fls{£), and thus {Ms{e))^ ^ (^'ff(^S))~^ Fs and {Ms{e))^^Hs have subsequential limits / and h = -Re J f"^, holomorphic and 
harmonic in ^ls{^)i respectively. An important observation, proven in Lemma 3.10 below, is that h cannot be identically 
zero. This implies that for any < ci < e we have Ms{ei) < CMg{e) with some constant C independent of 6. 

Hence, applying the diagonal procedure, we may assume that {Ms{e))~ ^ {'&{d))~^ Fg tends to a limit / (and {Ms{e))~^Hs 
tends to h = J f^) uniformly on each ^^(ei), < ei < e. Arguing as above, we see that h is harmonic in 
f2\{a, ai, . . . , afe}, satisfies Dirichlet boundary conditions, has positive outer normal derivative, and is bounded from 
below near ai, . . . , afc. Moreover, the last step of the proof given above shows that 

h(^^ = \imiMs{e))-\d{S))-^m f (Fsiz) - F[Cs .a]{z)f dz = limiMsie))-\^id))-^m ( {Fs{z)fdz = h 

(5-s-O J S^O J 

is also bounded in a neighborhood of a (one can neglect i^[(Cj.a](^) in the last expression since (i?((5))~^F[c^_(j](2;) tends to 
/[C,o] and Ms[e) — >■ oo). By the maximum principle, this yields h = 0, which is a contradiction. □ 

Lemma 3.10. In the notation of the proof above, none of the subsequential limits of {Ms{e))~^Hs is identically zero. 

Proof. Suppose by contradiction that {Ms{e))~^ uniformly on compact subsets of fi5(e). Let z^^^ be the point 
of ^^^(e) where the maximum Mg{e) of \Hs\ is attained. By sub- and super-harmonicity of H° and H* respectively, z™'^ 
belongs to one of the discrete "circles" ^o(e) := : e < — a| < e -|- 5^}, la^ := {z : e < \z — aj\ < € + 55} (recall that Hg 
vanishes on the boundary). By passing to a subsequence, we can assume that z™'^^ ^max_ 

Suppose that z"^'^^ G £aj (e) for some 1 < j < k. Since H* is superharmonic inside £aj (2e) and tends to zero uniformly on 
£aj (2e) by our assumption, its minimum m^(2e) inside iaj (2e) also tends to zero. Thus, z™'" g and Hg{z"^^^) = Ms{e). 
Hence, by subharmonicity of Hg, we can find a discrete nearest-neighbor path 7° := {0™^ = Zi ^ Z2 ^ ■ . ■} C Vq^ with 
Ms{e) < • • • < Hg{zi) < iJ|(zi_|_i) < . . . , which may only end up at Oj, where the subharmonicity fails. By Remark 3.10 
in [ChSml2], the functions H' — mg(2e) and — mg{2e) are uniformly comparable at adjacent points inside la^ (e). Since 
m^(2e) ^ as 5 ^ 0, while Ms{e) 00, we have H*{z) > cMs{e) for some absolute constant c > and all z G 7*, where 
7* is the set of vertices adjacent to 7°. Due to the maximum principle (see Remark 3.7), we infer that, for any z, 

Hl{z) > cMs{e)hm;{z) + (1 - hm;(^))m^^(2e), 

where hm*(z) is the discrete harmonic measure of 7 viewed form 2; in {2; : |z — a| < 2e}. But if z is close enough 
(independently of 6) to z"^'^^, then, by discrete Beurling estimate, hrRs{z) is close to one. This shows that, for such z, 
{Ms{e))~^H*{z) > c' > for an absolute constant c', and hence {Ms{e))^^Hs{z) cannot tend to zero in 0(e). 
It remains to treat the case z"^'^ G ^a(e). Consider the function 

(M,(e))-ifff ) = {Ms{e))-\mr^^^I{F5{z) - Fyc,,a^{z)?dz. 

Note that it tends to zero on £aj(2e), since one can neglect the term F^£^^a\{z) (recall that (i9((5))~^F[Cj,a](-2) — >■ /[c,a] and 
Ms{e) 00). Therefore, by Remark 3.8 and the maximum principle, it also tends to zero on ^aj (e), so we consequently 
derive that each of the functions 

Fs - F[C,,a] Fs Hs 

{Ms{e))yH{6)' Ms{e))yH{5) Ms{e) 

tends to zero uniformly on (aji^)- In particular, 1 = {Ms{e))~^\Hs{zf^'^^)\ — )• 0, which is a contradiction. □ 

3.5. Analysis near the singularities. We now pass to the most delicate part of our analysis: matching the second-order 
terms in the values Fs{a+ ^) with the second coefficient in the expansion of the continuous spinor near the "source" a. 

Proof of Theorem 2.17. Let TZa denote the reflection with respect to the axis {a + x,x gM.} and As be a small neigh- 
borhood of a in O^ n TZa (^5)- Denote by A^ the sheet of a double cover [A^, a"*] \ La such that ^^(a + ^) > on this 
sheet, where Fs = F[Q^^a;ai,...,ak] and La = {a + x,x < |} as before. We define a real-valued function Ss : ^ M as 

1 fFs + Fj^^ \ 
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where F^^^ = -F'[K„(n5),a,K<,(oi),...,K„(afc)]' = ^[o,a,ai,...,Ofc] and the full-plane functions -F[Q,a], were constructed 

in Section 3.2. Recall that F[Cs,a] («+¥)= 1' <^[Q,a] ^ + = ■^^('^ + ¥) by symmetry. Thus, 

Ss{a + f ) = (i? (<5))-i (^^(a + f ) - 1 - 25ie^ • (5) (3.19) 
and our goal is to estimate this value. Note that Sg vanishes on the cut both -Fp^.o] and G[Cs,a] vanish by construction, 

and = —Fg on La due to the spinor property [Fg changes the sign between opposite sides of La, since they belong to 
different sheets). Note that Ss is discrete harmonic in \ {^^ + ^} since all terms are discrete harmonic, and discrete 

harmonic at a + ^ also, due to Lemma 3.5. Therefore, for any (small) r > 0, discrete Beurling estimate (3.5) implies 

Ss{a+^)<C5^ ■r-^uyayi\Ss\, (3.20) 

where iair) denotes the "discrete circle" {z : r < 1^; — a| < r + 5(5}. It follows from Theorem 2.15 and Lemmas 2.13, 2.16 
that, for any r > 0, one has 

Ss s:=Jie[i(/ + /(^))-/[c,„]]-23?e^-5[c,a], 

uniformly near the circle {z ■.\z-a\= r}, where / = /[n,a;oi,...,afcl, /^'^^ = /[7^„(0),a;K„(al),...,K<,(afc), /[c,a](^) = l/s/z-a 
and g[c^a\{z) = — a. Recall that, by definition of the coefficient A, one has 



/ - /[c,a] = ^-Ay/z - a + 0{\z - a| 2 ), 2; 



a. 



It is easy to check that f^^\z) = f{TZa{z)) (since this spinor solves the corresponding boundary value problem), hence 

/^'^^-/[Ca] =2:4^/?^ + 0(|^-a|i), z ^ a. 
Thus, we arrive at s{z) = 0{\z — a\i) as z ^ a, which means 

max\Ss\ — > max \s\ = 0(r). (3-21) 

^air) 5->0 

Combining (3.19)-(3.21), one concludes that, for any given r > 0, 

Fs{a+ f) -l-2m A- 6 < C^{d)6ir, 

if 6 is small enough. Since '!?((5) = 0(6^) by (2.16), and r can be chosen arbitrary small, this yields (2.18). All estimates 
are uniform with respect to a, ai , . . . , afe at definite distance from the boundary and each other. □ 

To prove Theorem 2.19, we need to analyze the discrete spinor -F[[2,5,a:h] near the point b. In contrast to Theorem 2.17, 
where second-order information near a was extracted, here we only need to match the first-order coefficients. Note that the 
situation is slightly different from the first-order analysis near a, as F^Q^ a-b] does not have an explicit discrete singularity, 
remaining s-holomorphic near the branciiing b, while its limit blows up at b. Still, the strategy for the identification 
resembles the one used above and appeals to a symmetrization with respect to the horizontal line passing through b. 

Proof of Theorem 2.19. Let TZb denote the refiection with respect to the axis {b + x : x and A5 C fi^ fl (fi^) 
be a small neighborhood of b. Fix the sheet = A5 \ L;, on a double cover [O^, a; b] so that Fi{b+ ^) = iBs with Bs > 
(recall that the values of discrete spinors on V^^ are purely imaginary), and do the same for the continuous spinor: fix a 
sheet so that f[n^a;b]{b + x) > as x \.0 (see (2.12)). 

Let Ws := hm|j,_j_|j( • ) : \ — >^ M denote the harmonic measure of 6 -|- | in the slitted discrete plane \ Lf, and 

where Fg = J^[n^,a;b] and F^'^'' = i^[7^„(o),K6(a);b]• By symmetry, f]^^(6 -h |) = iBs (if one fixes the sheet for F^^^ by the 
same condition Fg''^\b -I- |) G iM+). Note that the function Ws is harmonic in A^ and vanishes everywhere on the cut 
Lb, including the point b + |: for 2; e Lb \ {6 + |}, one has Fs + F^^ = = Ws due to symmetry reasons and the spinor 

property of Fs, and Ts{b + ^) = since Ws{b+ ^) = 1. 

We have to prove that Bs ^ B = B[Q^a:h] as S ^ 0. By passing to a subsequence, we may assume that Bs ^ B for some 
B e [0, -|-oo]. If B is finite, then Theorem 2.15 and Lerrmia 3.3(iii) give 



Ts{z) — > t{z) :-- 

(5->0 



(3.22) 



uniformly on compact subsets of {2; : < - 6| < r}, where / = /[Q,a;6] and f^'^'^{z) = /[7Ji,(n),TC6(a);6](^) = -f{T^b{z)) on 



the sheet Aj". Note that 



t{z) = i{B-B)-m^^=+0{\z-b\i), z^b. (3.23) 
\/z — b 
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Therefore, B ^ B contradicts to the maximum principle for the discrete harmonic function T^: (3.22) and (3.23) imply 
that, for sufficiently small (5, its values inside of the circle ibif) around h are bigger then on (.b{r). 
The similar argument works, if S = +oo: in this case one would have 



«-!-o s/z-h 



which contradicts to the maximum principle for Ts, just as before. 



□ 



A. Appendix. Spinors in the half-plane and decorrelations identities 

In this Appendix we explicitly construct the continuous spinor which satisfies Definition 2.8 in the half-plane. Recall 
that, along with (2.6)-(2.8), one should also impose the condition /[H,o;ai,...,afc](2^) = 0(|2;|~^) at infinity, so that after 
a conformal mapping to a bounded domain with smooth boundary the spinor remains bounded. We argue that the 
construction of /[H,o;ai,...,ofc] boils down to a system of k+1 linear equations, which is always non-degenerate. Consequently, 
we derive the formulae for A\^^a;ai,...,ak] ^^^^ by integration, the continuous correlation functions. Although for small k 
it is possible to check by ugly brute force computations that the functions we get coincide with those predicted by the 
Conformal Field Theory arguments (recall that the computations of 1-point and 2-point functions were done in Section 2.7), 
we prove this coincidence only for some particular configurations when all a, ai, . . . , ak lie on the imaginary axis. This 
turns out to be sufficient to prove Lemma 2.22. 

A.l. The explicit formulae for spinors. Let us introduce some notation. For w e H, denote 

Pw{z) := {z-w){z-w). 
Given a,ai, . . . ,ak € H, we look for f[m,a;ai,...,ak] the following form: 

P{z) 



fp{z) := e 4 



\/Pa{z)PaA^)---PaA^y 



where P{z) = X]n=o 'Z"-^" ^ polynomial of degree k with real coefficients. Clearly, Jp is a holomorphic spinor on 
[H, a; ai, . . . , ak] which satisfies the boundary condition (2.6), hence we have to determine the coefficients qo,qi, . . . ,qk 
from the conditions (2.7) and (2.8). Since / \/Pa^{z) = (29mam)(-z — am) + 0{\z — am\^) as z ^ am, these 



conditions can be rewritten as 



J2 Mo^n{a; ai,..., ak)qn = 1, 

k 

Mm,n{a] Ol, . . . , ak)qn = 0, 

n=0 



(A.1) 



m 



1, . . . , fc. 



where 



< m, n < fc 



Mm,n = Mra,n{a;ai,...,ak) ■■= 3ffe 
(here and below we use the notation oq = a). If this linear system is non-degenerate, the Cramer's rule implies 

(In 



qn{a;ai,...,ak) = {-iy det [Mm,s\- D \ where D= det [M„ 

m^O^s^n ' 0<m,s<k 



(A.2) 



Proposition A.l. For any k + 1 different points a,ai, 
is given by (A.2). The spinor 



,ak € H, the linear system (A.l) has a unique solution which 



f[C,a;ai,...,ak]iz) ■= e * 



-Pa;ai,...,afc {z) 



\/Pa{z)PaAz) ■ ■ -Pa^zY 



-Pa;oi,...,Ofc (-2) '■— ^ ^ (jn^ i 



(A.3) 



n=0 



solves the boundary value problem (2.6) (2.8) and satisfies the additional regularity condition /[c,a;oi,...,ofc] (-z) = 0{\z\ ^) 
at infinity. Moreover, all the spinors /[c,a;ai,...,afc](-z) a'^e uniformly bounded, if z is e- away from a,ai, . . . ,ak, and 



AMia;ai,...,ak) 



1 



4(a - a) 



1 



a — Os 



a — as 



dzPa:ai,...,ak (z) 



2-Pa;ai,...,afc {z) 



(A.4) 



Proof. On the contrary, let the system (A.l) be degenerate. Then, a nontrivial solution to its homogeneous counterpart 
gives rise to a spinor fp that satisfies (2.6) and (2.7), as well as [lim^^a y/z — afp{z)] = 0, and thus vanishes identically, 
since 



< 



l\fpiz)fdz = 2TTlim{z-a){fp{z)f + 2ny] 1™ {z - as){fp{z)f < 0, 

Jm. ~l z^as 



(A.5) 
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which is a contradiction. Thus, (A.l) is non-degenerate and has a solution (A. 2). Let f[c,a;ai,...,ak] be given by (A. 3), thus 
satisfying the boundary conditions (2.6), and 

A = A(a;ai,---,afc) := hm^^.av'i^-a/[c,a;ai,...,ofe](^), < s < A;. (A.6) 

Note that the Unear system (A.l) means 5fte /?o = 1 and 5Re /3s = for all other s, hence the conditions (2.8) are satisfied. 

The contour integration argument (A. 5) ensures that /3q — (|/3ip -I h |/3/cP) > 0. In particular, /3o should be purely real 

(and so /3o = 1), i-e. the condition (2.8) is satisfied too. Moreover, one has 

h[C,a;ai,...,ak]iz) ■= / {f[C,a;ai,...,a„]{z)fdz = Gm{z\a) - ^0lG^{z\as), 

where Gm{z\ a) is the Green's function in the half-plane. Thus, (/[c,a;oi,...,oii](^))^ = \ [dzGm{z; a) - X^^^i PldzGm{z; as)] 
are uniformly bounded, since the gradients of the Green's functions are uniformly bounded, if z is e-away from the 
singularities. Finally, the formula (A. 4) readily follows from (A. 3). □ 

A. 2. Conformal Field Theory predictions and the particular case a,ai, . . . ,ak G It was predicted by the 

methods of Conformal Field Theory [BuGu93] that in the continuum limit, the spin correlation functions take the following 
form: 

k 



(Coo • • • faji )lHI,CFT ~ 



1 



=o(23ma™)« 



E n 









as 




2 


as 







(A.7) 



Remark A. 2. Recall that we have considered two simplest cases fc = and fc = 1 in Section 2.7. Namely, we found 
the explicit solutions (2.20) and (2.21) to the boundary value problem (2.6)-(2.8), computed their asymptotics near a 
and obtained the corresponding coefficients ^H(a) and An(a]ai). Then, one can directly check that those coefficients 
coincide with the logarithmic derivatives of the quantities (A.7), thus establishing this Conformal Field Theory prediction. 
Unfortunately, this procedure becomes much harder as k grows. Having in hand explicit formulae given in the previous 
Section, one can compute »4H(a; ai, . . . , a^) and integrate the answer, but it is not easy to check whether it coincides with 
(A.7) or not. Note that the case A; = 2 is still doable by hand. In this case, it is easy to check that, for a proper real 
constant C, the spinor 

, , . _ ^ in. C(r2Pa,{z) + (:oij;„,(,:) - f.l2j^„„(.--) _ _ 

/pH,ao;oi,a2] l^J ' / , . , . = ) Csm ■— \as — am\\as — ttml, 

satisfies the conditions (2.6) (2.8), which allows us to verify the CFT prediction (A.7) for the 3-point fimction by brute 
force computations. In the next lemma we prove (A.7) for all k in the special case when all point are purely imaginary, 
which is sufficient for our purposes. 

Lemma A.3. // = iwm e iR+ for all m = 0, . . . ,k, then d^^ log((Too • • • '^afe)H,cFT = ^-^H(ao; ai . . . afe). 
Proof. We can rewrite the function (tTao ■ • ■ ^^0^)11 cft given by (A.7) in the following form: 

2-|(fe+i) 



E 



n 



ami ^ la, -ami " 



Observe that the logarithmic derivative of the first factor exactly matches the first two terms of (A. 4). Hence it suffices 
to show that 



dwo log 



|Ws - Wml ^ \Ws +Wml = 



idzlogPao;ai,...,adz)\, 



- jUo=±l,...,([ife=±l s<m 

= dn,\ogPao;a„...,ad^w)\^^^^. 

Let, say, wq < wi < ■ ■ ■ < Wk, the other cases are treated by similar computations. Then, 

E \Ws-Wm\ ^"2"" |Ws + "2""' = ^ (Wm - MsMm^s) 

jUo=ilv)/^fc==tl s<m /:Ao==tl5-",/^fc=il s<m 

k k 

E WlJ'Z^il^n.Wm- l^sWs) = ^ [J /X^^J ^^dct^^ [(/Xr»W„)"] 

jUo==tl v)/^fc==tl m=0 s<m /Lto==tl,...,;y.fc=±l m=0 ~ ' ~ 

y dct = det [(l + (-l)™+")u;;j^]. 

Aio = ±l,...,A'ic=±l 

On the other hand, Pai{am) = sign(/ — m) • \wi — WmWwi + Wm\, thus 

Mm,s = [{iwmy ■ i'^Rm] , where Rm = 11;^™!"^; - Wm\\wi + w, 



(A.8) 



(A.9) 
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are real and do not depend on n. Note that SRe = 5(1 + Therefore, (A. 2) reads 



2 ^ 



g„ = D-i.(-l)" det [M„,,] = C •(-!)" det [1^+^(1 + (-ir" ■ ^ 



= C- det [(1 + 

= (-l)^C.Hr det [(l + (-l)™+>^], 

where C = D-'^2-^Ri ...Rm, and 

fe 

Pao;a„...,a.M = ±C • ^ det [(1 + (-1)™+'X] = i^C • ^ det [(1 + , (A.IO) 

n=0 

since detm^o,s#n [(1 + {—^)"^"^)xm,n[ = for odd n and any aJm^^. Thus, (A. 8) readily follows from (A. 9) and (A.IO). □ 



A. 3. Decorrelation identities and proof of Lemma 2.22. We start by proving that the half-plane spinors /H,a;oi,...,Ofc 
behave in a continuous fashion, when one of the points Oi, . . . , a/j approaches the real line and "disappears" there. 

Lemma A.4. As dE., one has f[m,a;ai,...,ak\{^) /[ IH[,a;ai,...,a/E_i] 

{z) andA[Q, 

,a;ai,...,afc] ~^ -^[n,a;ai,...,ofc_i] > Uniformly 

with respect to the positions of points a,ai,..., ak-i,z, provided they are at least e- away from the boundary and each other. 

Proof. This can be easily proven by compactness arguments, but we prefer to give an explicit construction which also 
allows one to estimate the convergence rate, if necessarily. Denote 

9 ■■= /[H,a;ai,...,afe_i] {z) ' {z), where Ta^ (z) := -j===== 

\/ {z - ah){z - ak) 

By definition, / satisfies the spinor property on [C, a; ai, . . . , Ofe], the boundary conditions (2.6) and g{z) = 0{\z\~^) at 
infinity, but the conditions (2.7), (2.8) fail since the values riy{as) are nonreal. Let = /3s(a; ai, . . . , ak-i) be defined 

by (A. 6), and 



/3o 



lim^^a ^z-ag{z) = 15^^'^^ ■ ra^{a) = ra^ia), 
lim^^a, \/z - as g{z) = pi'''^^ ■ raja^) 
lim;,^Q, y/z - at,.g{z) = /[H,a;ai,...,afe_i](afe) • (iSjmafc)5, 



where 1 < s < — 1. Note that g can be represented as 

k 

9 = 5Z ^^'^s ■ /[H,a,;ao,...,a,,...,afc]- 

Indeed, the difference of two sides would solve the homogeneous version of the boundary value problem (2.6)-(2.8), which 
is impossible since the solution of this problem is unique due to the contour integration argument (A. 5). Therefore, 

fc 

9- /[H,a;oi,...,afe] = (^^^ /^O - l)/[H,a;ai,...,afc] + ^ • /[H,a,;ao,. ,afe] • 

s=l 

Note that /3o = ra^ia) — >■ 1, 3?e/3s = — 9tn/3s'^~^''9mra^(as) — and /3fc — )■ as y — > 0, and all the spinors in the 
right-hand side are uniformly bounded due to Proposition A.l. Thus, we conclude that 

f\a,a;ai,...,ak]{^) " f[U,a;ai,...,ak-i]{z) = [f[M,a;ai,...,ak]{z) " 9{z)] + [1 " rak{z)]f\a,a;ai,...,ak-i]{^) ^ 

since both terms vanish as y — )• 0. The convergence of A[Q^a;ai,...,ak] follows by the Cauchy integral formula. □ 
We are in the position to prove Lemma 2.22. 

Proof of Lemma 2.22. We first define the correlation functions {aaCTa^ ■ ■ ■ (^a^)^ ™ half-plane. For A; = 0, 1 we use 

the explicit formulae (1.3), (1.4), and we define them to be cxp[J' £n,fe+i] in the general case, where the differential form 
Cii^k+i is given by (2.23). Note that this form is symmetric with respect to a, ai, . . . , Ofe and exact on the manifold f2'^+^ 
of all [k + l)-tuples of distinct points in O due to Proposition 2.20. Therefore, the function {<Ja<Jai ■ ■ ■ <^ak)m symmetric 
with respect to a,a-i, . . . ,ak and well-defined on fl'^~^^ up to a multiplicative constant which we fix so that 

{(TaCTai ■ ■ ■ (Tak)u = i'^a.cra.i ■ ■ ■ (^ak)m,CFT a, tti , . . . , ttfe € «+. (A.H) 

As it was established in Lemma A. 3, this normalization (everywhere on the imaginary line) docs not contradict to our 
definition of the correlation function as an exponent of the primitive of Cn^k+i- Therefore, (A. 11) holds true for any 
a,ai, . . . ,ak lying on the same vertical line: indeed, the multiplicative normalization cannot depend on the horizontal 
position of this line, since for such configurations the coefficients Auiao', ai • • • ctfc) are purely imaginary due to Lemma A. 3. 
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Let a — )• 9EI along some vertical line, ai , . 
say, hs := Jieo + is. Write 

{'^a^ai ■ ■ ■ f afc)H 



, afe e H and 61, . 



, hk be some fixed points on the same vertical line as a, 



'ifc/H 



Since a, 61, . . . , 6^ are on the same vertical line, we know from the explicit formulae (A. 7) that the first term tends to one 
as a ^ 0. Also, we have 



log 



+ 1 

H 



\<7a<7bi 



(a,ai,...,ak) 



{ai,...,ak) 



-a,fe+l 



{a,bi,...,bk) {bi,...,bk) 

According to definition (2.23), the differential forms Cn_k+i and C^^k consist of the terms jRe [An{zj; zi, .... Zj, .... Zk, a)dzj] 
and jjfe [Anizj; zi, .... zj, Zk)dzj], respectively, which are uniformly close to each other due to Lemma A. 4, as (zi, . . . , Zk) 
runs from (ai, . . . , a^) to (61, . . . ,bk) and ?sma 0. Note that, since the convergence is uniform, it does not matter, 
whether a dW along a vertical line or not. 

Thus, we have defined the half-plane functions {(JaCTai ■ ■ ■0'ak)m ^^^^ decorrelation identities (2.27) are fulfilled. We 
then define the correlation functions {<Ta<Tai ■ ■ ■ ^ak)ri ^ given bounded domain Q using a conformal map ip -.0. and 
the covariance rule (1.2). The computation given in Remark 2.21 shows that, being defined in this way, \og{(Ja(yai ■ ■ ■ '7ak)n 
is indeed a primitive of the form CQ„k+i- Since the covariance rule (1.2) implies the conformal invariance of the ratios 
(o'aO'ai • • • <7ak)Q/{(^a)Q{<7ai ■ ■ ■ o'afc)Q, One obtains dccorrelation identities (2.27) as a — )• dfl (which means ip{a) — )• H). The 
second part of (2.27) was proven in Remark 2.23. 

It is easy to see that conditions (2.27) fix the unknown multiplicative normalization of the correlation functions uniquely: 
the full-plane normalization fixes the two-point function and decorellation identities inductively fix all others. Therefore, 
the functions {<Ta<Tai ■ ■ ■ (7ak)n constructed above do not depend on a particular choice of the conformal mapping ip. □ 
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